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ABSTRACT 

Liquid viscosity has been included in the Bellman-Pennington 
theory of the Taylor wave in a liquid vapor interface. Predictions 
of the "most susceptible" wavelength, and of the wave frequency, 
are made as a function of a liquid viscosity parameter and the Bond 
number. 

The stability of a gas jet in a viscous liquid has been studied 
and the ensuing result is used to predict the peak heat flux on 
large horizontal heaters. 

Experimental measurements of the dominant Taylor wave and 
its growth rate are made during the film boiling of cyclohexanol on 
cylindrical heaters. The results bear out the predictions quite well. 
The thickness of the vapor blanket surrounding a cylindrical heater 
has been measured and a correlation suggested. The effect of large 
fluxes of vapor volume on the dominant wavelength has also been 
noted . 

The peak heat flux is observed on relatively large flat plate 
heaters both at earth normal gravity and at higher gravities for 
boiled liquids that are nearly inviscid ( e.g . acetone, benzene, metha- 
nol isopropanol and water) and for boiled liquids that are significantly 
viscous (e.g . cyclohexanol). The peak heat flux is also measured 
when cyclohexanol is boiling on cylindrical heaters. Theoretical 


iii 



results of the peak heat flux are compared with the experimental 


data. 

The effect of finite geometry of flat plate heaters on the peak 
heat flux is also discussed. 
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I. INTRODUCTION 
Brief History of Previous Studies. 

Nukiyama [I] 1 , in 1934, otserved the stable regions of the 
now well-known heat flux versus temperature difference curve for 
pool boiling. He identified these regions as nucleate boiling and 
film boiling respectively. Nukiyama was also able to locate the 
two extremes of this curve, i.e. the peak and minimum boiling 
heat fluxes. Drew and Muller [2] subsequently observed the in- 
termediate region which is now called the region of transitional 
boiling. 

Since then, extensive studies have been made of the various 
regimes of Nukiyama’s boiling curve. These studies have been 
mainly motivated by the association of the boiling process with 
much larger heat transfer coefficients than those occurring in 
other processes. In this chapter we will briefly review the evo- 
lution of the hydrodynamic theory as applied to the peak and min- 
imum heat flux on flat plates and cylindrical heaters. 

In the late forties, Kutateladze [3] suggested that the transi- 
tion from nucleate to film boiling was purely a hydrodynamic 

1 Numbers in square brackets represent entries in Reference section. 
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phenomenon. Based on similarity considerations related to the 

equations of motion and energy, he identified two dimensionless 
2 

groups : 


P 


q 

max 


h U 

g fg g 


and 


°(P f - P g )9 

4 2 

U p * 

9 9 


(D 


where q is the peak heat flux and U is the velocity of the 
max g 

gas in vapor jets. Combining these two groups he gave a very 
successful correlation for the peak heat flux: 

fq ^ • . . * o.i6 nr h IW-p )T /4 

v. maxxinviscid v 9 f9L ^ f g J (2 


Kutateladze did not suggest any restrictions as to the heater geo- 
metries that equation (2) should describe. 

Borishanski [4] correlated peak heat flux data obtained by 

various authors for circular heaters and horizontal strips. In a 

3 

personal communication , Borishanski recalled that the data he 
used in forming the correlation were generally free from induced 
convection effects. His modified equation, similar to equation (2), 
with the effect of liquid viscosity included, was: 

Vax F = '°' 13 + ^ h fg[°«Pf-0g>]' /4 } < 3 > 

2 

Symbols commonly used are explained in the section on nomen- 
clature. 

3 

Conversation with J.H. Lienhard in Minsk, U.S.S.R, May 17, 1972. 
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where 


N s 





(4) 


Costello et al. [5J showed that the peak heat flux could be 


strongly affected by induced convection. Subsequently, Lien hard 
and Keeling [61 correlated their peak heat flux data for ribbons 
subjected to induced flow from the sides on a i/nlJ axis. In 
their case, was the characteristic width of the ribbon: 


L’ = «-^g(p f -p g )/a 


(5) 


In 1958, Zuber and Tribus [7] explained that the transition 
from nucleate to film boiling was a consequence of hydrodynamic 
instability. The peak heat flux occurred when vapor jets became 
Helmholtz unstable and hence collapsed; thereby obstructing the 
flow of liquid toward the heater. By making certain assumptions 
regarding the spacing of the jets, their diameter, and the wave- 
length at which the vapor jets became unstable, Zuber obtained 
an expression for the peak heat flux from a horizontal surface: 

tt r -1/4 

q ma* 2 " Sl/% V 0 *V p g>] (6:i 

Equation (6) applies as long as the system pressure is not close 
to the critical pressure. We have not amplified the derivation of 
equation (6) here as its development will form a part of our dis- 
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cussion in Chapter IV. 

Cichelli and Bonilla [ 8 ], Berenson [9], and Costello et al. 
[5] measured the peak heat flux on flat plate heaters in the ab- 
sence of any induced flow from the sides. The limited data of 

the first two authors for clean surfaces seem to bear out equa- 

4 TT 

tion ( 6 ) with a little change in the constant , ^4 • However, 
Costello's data for distilled water are much lower than predicted 
by equation ( 6 ). Thus, more data must be obtained for flat plate 
heaters with different liquids and at various gravities before one 
can judge whether or not equation ( 6 ) really represents the peak 
heat flux on horizontal surfaces. Part of the present effort will 
be aimed at doing this. 

In 1964 Bobrovich, Gogonin and Kutateladze [ 10 ] and in 1965 
Lienhard and Watanabe [11] independently showed that diverse 
peak heat flux data on cylinders of various radii and at different 
gravities could be correlated by: 


q = q • f(R') 
max max.. 


where f(R') is an undetermined function of R' and 


(7) 


R' = R^/g(p f -p )/o 

y 


( 8 ) 


Equation ( 8 ) may be taken as a dimensionless parameter combin- 


4 

This is explained in detail in Chapter IV 
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ing the size of the heater and the gravity acting on the system. 

Later, in 1970, Lienhard and Sun [1 2] proposed a hydro- 
dynamic theory of boiling for cylindrical heaters and gave the 
following equation which was verified by about 900 c 1 rnax data over 
a wide range of heater size, liquids, pressures, and gravities: 

^max^inviscid ^max^ *£0,89+ 2.27exp(— 3.44-v/K*)] 

1 /4 

for R'^0.15. This study conclusively proved that q rnax '“9 fcr 
R'^3.5 for all gravities. However, for small R', the peak heat 
flux depended on gravity not only as it appeared in q max but 
also through the functional dependence upon the dimensionless 
radius, R'. Thus the disagreement arising out of Costello and 
Adam's [13] observation that the peak heat flux did not depend on 
one-fourth power of gravity at lower gravities was resolved. 

Zuber also realized that the steady cyclic release of bubbles 
at discrete locations in film boiling was nothing but a classical 
example of Taylor Instability. He argued that at the minimum 
heat flux, the spacing or the bubble releasing nodes lie between 
'most susceptible" and "critical" Taylor wavelengths and that the 
average rate of penetration of vapor region into liquid was slower 
than that predicted by Taylor instability. Zuber assumed that 
the bubble diameter at the departure was half of the dominant 
wavelength and that the bubble grew linearly in accordance with 
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Taylor instability theory up to a height equal to 40 percent of 
wavelength. He based his arguments on the experimental obser- 
vations of Lewis [15] and Allred and Blount [16]. Zuber aver- 
aged the growth rate during linear growth and for the bubble 
frequency obtained: 


where 


f 1 dT ^ 

b X . dt 

d 


— 0.4X 

dig _ 1 f* _dig 

dt “ 0.4X J dt 
o 


dig = 0.2 id 1 
d d 


for X = 


X 


d 


Thus, for the minimum heat flux on flat surfaces he obtained: 


(q_ 




( 10 ) 


Berenson [9] subsequently made several observations of the 
minimum heat flux on flat plates. He argued that Zuber should 
have time-averaged the growth rate rather than averaging over ig. 
A considerable portion of wave growth may be nonlinear, hence a 
complete growth history is needed to form a proper time average. 
Berenson himself could not form the average because of the 
limited information available to him. However, he suggested 
that the experimental constant in equation (10) should be 0.09 in- 
stead of 0.177. 

In 1963, Lienhard and Wong [17] incorporated the effect of 
transverse curvature of the interface in the instability theory. 
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They were thus successful in predicting "most susceptible" Taylor 
wavelength to be expected in film boiling on horizontal cylinders. 
Based on arguments similar to Zuber’s they obtained the follow- 
ing expression for the minimum heat flux: 


( q min) inviscid 


C q minp)j nv j SC | C | 


r 1 ♦ 289 

’[ 2 2 
R' (2R» +1) 



(ID 


where d m inp * s 9^ ven by equation (10) with the constant changed 
to Berer.son's value. Their experimental observation for alcohols 
and the data of Kovalev [18] for water confirmed the trend given 
by equation (11), but the data of the two authors define different 
constants in equation (11). The difference was due to the way 
the wire was mounted during the experimental observations. The 
end mountings can have a significant cooling effect on the wire 
and may cause an early collapse of film boiling at the ends. 

In 1967, Baumeister and Hamill [19] made a theoretical 
prediction of the heat transfer coefficient for film boiling from 
horizontal wires. They solved the equations of motion and energy 
and optimized the heat transfer from the wire. Their rough model 
of the film boiling process also enabled them to predict the vapor 
blanket thickness around the cylindrical heater. 

Later, in 1969, Lienhard and Sun [20] measured wavelengths 
during film boiling on horizontal cylinders. Their experiments 
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included a wide range of liquids, wire sizes and gravities. For 
a particular dimensionless radius, R’, their data for wavelengths 
showed a significant variability. However, it was found that all 
the wavelengths were consistent with frequencies within ± 10 per- 
cent of the maximum frequency. 

The present study will concentrate both on the "infinite" flat 
plate and cylindrical heaters. The induced convection situation 
will be avoided in the first case by limiting the configuration of 
the boiling vessel to that of a flat plate with vertical side walls 
to approximate an infinite flat plate. In the second case, the 
geometry itself is such that the induced flow is free to pass by 
without interfering with the bubble flow path. 

Present Objectives. 

Nearly all the work done so far on boiling heat transfer 
ignores one of the important liquid transport properties, i.e. 
viscosity. This assumption may be valid at pressures that are 
not extremely low for almost all the Newtonian liquids usually 
used in boiling. Still, one would like to know the contribution 
of viscosity, if any. In this context we have the following goals 


in mind: 
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1) Taylor Wave Behavior , 

(a) Extend Bellman and Pennington's [14] Theory to deter- 
mine the effect of liquid viscosity on the ’’most susceptible" 
wavelength and che corresponding frequency. 

(b) Make experimental observations of the dominant unstable 
wavelength in film boiling on cylindrical heaters to pick up 
the effect of liquid viscosity, and make observations of the 
rate of growth of the unstable disturbance. 

Since, for the commonly available liquids, viscosity de- 
creases exponentially with temperature, the observations in 
part (b) will be made at very low pressures. Thus we 
would like to know: 

(c) The effect, if any, of large volume fluxes in stretching 
the dominant unstable wavelength. 

2) He at Fluxes. 

(a) Use above obtained expressions for Taylor wavelength 
and the wave frequency to discuss the minimum heat flux 
on an "infinite" flat plate and a cylindrical heater. 

(b) Attempt to predict how the peak heat flux varies with 
viscosity in both these geometries. 

(c) Put Zuber's equation (6) to a rigorous test by measui — 
ing the peak heat flux on an "infinite" flat plate heater at 
various gravities with different liquids being boiled. If 
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equation (6) does not represent the data for peak heat flux 
on an "infinite" flat plate, suggest a new theoretical model 
consistent with the observations. 

(d) Observe the peak heat flux on an "infinite" flat plate 
and a cylindrical heater, when the boiled liquids are viscous. 
Relate the data to item 2b. 



II. DERIVATION OF EXPRESSIONS FOR WAVELENGTH 
AND MINIMUM POOL BOILING HEAT FLUX IN 
VISCOUS LIQUIDS 


In 1950, G.I. Taylor [22] discussed the instability of the 
horizontal interface between two ideal incompressible fluids of in- 
finite depth. He showed that irregularities at the interface 
tended to grow if the acceleration was directed from the less 
dense to the more dense medium. Bellman and Pennington [14] 
extended his problem and showed how to take into account the 
interfacial surface tension and fluid viscosity. They gave closed 
form expressions only for the "most susceptible" wavelength 
( i.e. , the wavelength for which the growth rate for a small dis- 
turbance at the interface is maximum) and the corresponding fre- 
quency. They obtained solutions for the situation in which only 
the surface tension at the interface was considered, but they were 
unable to do so in the more general case in which both the li- 
quid viscosity and interfacial surface tension were considered. 

We shall begin with formulation of the instability problem 
similar to that made by Bellman and Pennington; we shall then 
incorporate the pressure contribution of curvature of a cylindrical 
surface as was done by Lienhard and Wong [17] for the inviscid 


11 - 
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case. Thus we shall attempt to obtain numerical values of the 
"most susceptible" frequency and corresponding wavelength as a 
function of liquid viscosity, as well as radius of curvature of the 
cylindrical heater. This evaluation in turn will be used to predict 
the minimum heat flux. 

Hydrodynamic Analysis. 

Before dwelling on the instability theory it is appropriate 
to state the various assumptions we are going to make: 

1) The fluids are incompressible. 

2) The fluids are Newtonian. 

3) The fluid depths are large as compared to wavelength. 
This is not exactly valid for the vapor blanket thickness in 
film boiling, but we shall show later that the above assump- 
tion is not unrealistic. 

4) The liquid is above a gaseous phase, such that the ac- 
celeration at the interface is always directed from the 
lighter to the heavier fluid. We assume this is so, because 
this is the situation encountered in all film boiling cases we 
will treat. 

5) Nonlinear effects are negligible. 

Figure la shows the interface between two incompressible 
viscous fluids of infinite depth. The linearized equations govern- 
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ng. la Interface between two incompressible 
viscous fluids of infinite depth. 



Fig. lb Schematic diagram of a typical 
configuration of film boiling on 
cylinders. 
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ing the motion are: 


u + v = 0 

x y 


u = - — p + — (u -*-u ) 

t p x p XX yy 


1 pi , 

v = p - g + — (v + v ) 

t p y p xx yy 


These equations are satisfied by 


and 


provided that 


and 


u = -to - \|i 
x y 


V = -0 + t 

y x 


p = p - pgy + p<p 
o t 


0 + 0 =0 
XX yy 


-(^ + t ) = 

p XX yy 


where u and v are the velocities in the x and 
spectively; p is the static pressure; p Q is the 
p is the density of the fluid. 

Let us take, for the liquid phase: 


A -Py+uut 

to = A e cos kx 

r f 


t = B e rr) f^' +a,t S i n k x 
f f 


( 12 ) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

y directions, re- 
total pressure, and 

( 20 ) 
( 21 ) 
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Pf = p o - P f gy + p 


( 22 ) 


where k is the wave number, a is the frequency, and: 


2 2 

m f = k + p f w/P f 


Similarly for the gaseous phase: 


ky+at 

(0 = A e cos kx 

g g 

m q y+at 

t = B e y sin kx 

g g 


and 


p = p„ “ p^gy + p^co ) f 
g o g g g t 


2 2 

m = k + P w/u 

g g g 


(23) 


(24) 

(25) 

(26) 


(27) 


The real part of m f or m^ has to be positive so the velocity re- 
mains finite as y becomes large. 

Considering that waves of height y = ri(x,t) are propagated, 
we obtain for the linearized kinematic condition at the interface: 


or 


= V f 


at 

r) = i) e cos kx 


where 

k(A f + B f ) 

i) “ 

o a 


(28) 


(29) 


( 30 ) 


The boundary conditions at the interface, i.e. , at y=r|, are: 



- 16 - 



- p f + 2 U f C Vy * -Pg + 2 w g ( Vy' < ” 1 ** 

MfCCVx +(u f J y ] ’ V C Vx + '"gV 


(31) 

(32) 

(33) 

(34) 


The substitution of equations (20) through (22); (24) through 
(26) and (29) in equations (31) to (34) gives four linear and homo- 
geneous equations in the four constants, A e , B_, A , and B : 

f f 9 9 


!<A. + m B — kA + m B =0 
f f f g g g 


(35) 


A + B + A - B =0 
f f 9 9 


(36) 


p9( p f - p g) k ok 


g(p--o )k ak" 


L ju 


U) 


2i a' •» 

i -2uk A + s 2|i km B. 

f fJfL uu ou ffjf 


+ [p u> + 2pt k ] A - [2pi km jB =0 
9 9 9 9 9 9 


(37) 


and 


A + w /k 2 + m 2 ) B + 2pi k 2 A - U (k 2 + m 2 ) B = 0 f38) 

rrr ft Q a Q qq 


The above equations have a non-trivial solution if and only if 
the determinant of the coefficient matrix is zero, i.e.: 




the "most susceptible" wavelength, while the effect on the 
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decrease it by a small amount. 

As proposed earlier by Lienhard and Wong [17], a two 
dimensional model can be used for cylindrical geometries. The 
effect of transverse curvature is treated in the form of an addi- 
tional oscillating pressure difference component across the inter- 
face, as shown in Figure 2. 

The expression for the transverse pressure may be written 


as: 


AP tr 2R 2 (1+d g /R) 2 11 


(40) 


where d g is the depth of the vapor blanket over the wire. 



Figure 2. Contribution of transverse curvature on a two 
dimensional wave. 



If we define a corrected radius, R * R[1+d /R], we can write equation (33) as: 
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[-g(P f -P )k + crk 3 +(p f +p g )u)J £p^k + m f )+p g (k + m g )J 

9 2R — 

c 

+ 4'JJkQi k + p m ] [p k + am J = 0 (39b) 

f g g g f r 

Our next step is to put equation (39b) in a more usable 
form, so that it can be explicitly solved for frequency and wave 
number, or wavelength. The vapor viscosity p^ is much less 
than p^. so it may be neglected in comparison to liquid viscosity. 
Thus we may write, for the frequency, x. 


«» = fiCPf+Pg* P f _P g» a ’ k * g> R C } ( 41 ) 

This expression relates eight quantities which are expressible in 
three dimensions. Using the Buckingham Pi-Theorem, we can 
recast '.his problem in terms of five dimensionless groups. For 
these groups we choose: 

(i) a dimensionless frequency 


1/4 

“■■[iT-] 

9 (P f -P g ) 


(42) 
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or a dimensionless wavelength. A, defined as 


A = 1/(>/3 K) 


(43b) 


(iii) a dimensionless liquid viscosity parameter 


M 



(44) 


The square of this group is very nearly the 
Borishanski number, N, as defined in equation (4). 


(iv) a dimensionless density 



P f +P g 


(v) a non-dimensional cylinder radius 


(45) 


-9( P f- P g) 1/2 


This number is related to the Bond number. Bo, 

2 

by Bo = R' . 

c 

Using these dimensionless numbers, we may write equa- 


tion (39b) as 
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(K 2 +fiM) 1/2 
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(K 2 +0M) 1/2 


ft 


r(K 2 +nM ) 1/2 


4ft K . K 

+ 2 1/2 

Mr 2Bo(K +QM) ' 


(47) 


Equation (47) gives a relation between frequency and wave 
number for an unstable disturbance (corresponding to Q positive 
and real), when the interfacial surface tension, transverse curva- 
ture and liquid viscosity are taken into account. 

From equation (47) it is clear that when the wave propaga- 
tion velocity is zero, there is no effect of liquid viscosity on the 
critical wavelength. Also, when M -• °® (i.e. , the liquid is invis- 
cid) equation (47) reduces to the same equation as obtained by 
Lienhard and Wong [16] for the inviscid case. 

We are interested in the ’’most susceptible" frequency, or 

the frequency for which growth rate of the disturbance is maxi- 

d £"2 

mum (i.e., the one for which — = 0). Differentiating equation 
dk 

dn 

(47) with respect to k and setting — — = 0, gives 

UK 
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4Q 


2 3/2 

(K +QM) 


r(K 2 +(7M) 3/2 MF 2 Bo(K 2 +QM) 1/2 


K 


2 3/2 

2Bo(K +QM) 


= 0 


(48) 


Thus we have two nonlinear equations (47) and (48) in two un- 
knowns, 0 and K. The two equations were solved numerically 
on an IBM/360 computer using a subroutine: XFNLES, which 
solves a system of nonlinear equations. This subroutine is avail- 
able in the Numerical Analysis Library [25] for the S/360, 
r was taken to be 0.9995 (^1) in the calculations. 

The "most susceptible" frequency and corresponding wave- 
length are plotted in Figures 3 and 4, respectively, as a function 
of the viscosity parameter M and Bond number. Bo. As evident 
from these figures, the effect of liquid viscosity is to increase 
the wavelength and to decrease the corresponding frequency. Fig- 
ures 5 and 6 show the effect of Bond number and liquid viscosity 
separately on the frequency wavelength relationship equation (47). 
Transverse curvature of the cylindrical heater reduces the "most 
susceptible" wavelength and increases the frequency. Figure 6 


shows that apart from increasing the wavelength, viscosity also 



- 24 - 



^J^oududnbej^diqjideosns 4 soui |( sseiuoisueuiip 
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dimensionless liquid viscosity parameter, M 
Fig. 3 Variation of £ 1 ^ with M for various Bond numbers 
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tends to increase the region of near neutral stabilitv slightly. 

By "region of near neutral stability" we mean the range of wave- 
lengths that can exist within, say, 10 percent of the maximum 
frequency. 


Minimum Heat Flux on a Flat Plate. 


We can make a prediction of the minimum heat flux by 
making an energy balance at the surface of the heater. 


q 

min r 


r latent heat transport" 1 

r number of bubbles! 

L unit volume 

unit time J 


[ 


volume of a bubble 


][ 


1 


area of the heater 


:] (49) 


or 


X 3 


q 

mirv 




(50) 


where we have assumed that the dominant wavelength is the 
"most susceptible" one, and the departing bubble diameter is 
half of it. Lienhard and Wong have shown experimentally that 
this is true within a certain range. In the above expression for 
q m inp» tbe bubble release frequency, f^, has been assumed to 
be buu , where b is an experimental constant. Here again we do 
not yet include the effect of viscosity on b. Experimental obser- 
vations will subsequently shed more light on the validity of these 


assumptions. Equation (50) may be written as 
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1/4 


(51) 


Dividing equation (51) by equation (10) we get 

q min 

Q . = = 1 .608 QA/T 7 

min (Pminc-)' . . . d d 

F ' min F mviscid 


(52) 


Minimum Heat Flux on a Cylindrical Heater. 


Following arguments similar to those for tne flat plate we 
may write an expression for the minimum heat flux on cylinders: 


4 V 

q rmn * C VW C W/ 2nR V 


(53) 


This may be rewritten as 

pjT 


g . =iJA 0 A 2 r— 2 1 

rmn 48 R d d Lg(p-p )J 

<■ g 


3/4 


1/2 


(54) 


Lienhard and Wong gave an expression for the inviscid 
case as 


(q . ). . 


. i p K 2 9(p c -p ) 
_ *> i_ g fg r___f_g_ 


g y/ 2 

''min'inviscid 48 -.3/4 R L p +p + (p +p )R‘ =J J 
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( 55 ) 


Dividing equation (54) by equation (55) we get 
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min 


'min (q . . . 

rruninviscid 


= 1.585 


Vd 

- 1/2 

i 


[1+1/2 3o] 


1/4 


(56) 


Thus, knowing Cl and A as functions of the liquid viscosity 
d d 

parameter, M, and the Bond number. Bo, one can predict the 
minimum pool boiling heat flux. As M is decreased, and 
vary such that the dimensionless heat flux, q . /(q . ). 

We shall return to question this particular result on pages 70-71. 


Conclusions. 

1 . The effect of liquid viscosity is to decrease the "most 
susceptible" frequency and to increase the corresponding wave- 
length. 

2. The effect of transverse curvature of the interface is 
to shorten the "most susceptible" wavelength and to increase the 
corresponding frequency. 

3. There is no effect of liquid viscosity on critical wave- 
length. 

4. The region of neutral stability (i.e. , the region close 
to the 'Ynost susceptible" frequency) increases slightly with liquid 


viscosity. 



III. EXPERIMENTAL OBSERVATIONS OF THE DOMINANT 
UNSTABLE WAVELENGTH AND ITS GROWTH RATE 
DURING THE FILM BOILING OF 
VISCOUS LIQUIDS 

An experimental program was carried out to observe the 
wavelength, its rate of growth, and the thickness of the vapor 
blanket surrounding the wire heater during film boiling in viscous 
liquids. Since the viscosity of most of the liquids is fairly low 
when they boil at normal pressures, the experime.its had to be 
performed at very low pressures. With the limitations of attain- 
ing maximum vacuum in the existing boiling apparatus in mind, 
a survey of various available chemtc •* ade. It was found 

that Cyclohexanol CH^CH^)^ CHOOH anc. ; ene Glycol 

CH (CHOH) CH OH were best suited for the However, 

nearly all the experiments reported here were performed with 
cyclohexanol as the boiling liquid. Various physical properties 
of cyclohexanol are plotted in Appendix B. 

A. Description of the Experiment. 

1 . Experimental Apparatus. 

A test capsule 8.9 cm wide, 8.9 cm high and 17.8 cm long, 
with glass windows in the 8.9 by 17.8 cm walls, was used in all 
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of the experiments. The capsule was made of brass and insul- 
ated on the sides with 0.63 cm thick styrofoam sheet. An elec- 
tric preheater and 0.63 cm diameter brass holders to support the 
test heaters were fitted to the capsule. Nichrome wires were 
used as test heaters and copper hooks were attached to their 
ends so they could be mounted in the brass holders. The hooks 
were attached in such a manner so as to minimize the end ef- 
fects. A 2.54 cm marker was mounted on the bottom of the cap- 
sule to provide a reference dimension for the reduction of photo- 
graphic data. 

Figure 7 shows a schematic diagram of the apparatus . AC 
power was employed in most of the experiments. The power 
supply to the wire was calculated from the measured current and 
voltage in the wire. A mercury manometer or a vacuum gage 
was used to note pressure inside the capsule. This pressure 
was corrected to take into account the head of liquid above the 
wire. An identical apparatus was used by Sun and Lien hard to 
make similar measurements. Full details are given in reference 
[ 20 ]. 

2. Procedure. 

A nichrome wire approximately 10 cm long was cleaned 
with soap and hot water to remove any grease or oily matter and 
then was rinsed with the test liquid. The wire surface was 
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Schematic diagram of the apparatus. 
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smooth anri had a cold-rolled finish (as supplied by the manufac- 
turer). The wire was mounted in the brass holders and the 
glass windows were positioned. The capsule was filled with 
reagent grade test liquid to a level of about 2.5 cm above the 
wire. The vacuum pump was started and the preheater was 
turned on to heat the liquid to saturation temperature. The pre- 
heater was turned off before energizing the wire. This was done 
to avoid any effects of convective currents and electric fields. 

The current in the wire was steadily increased until the 
peak heat flux was reached and the transition from nucleate to 
film boiling was observed. Thereafter the current was reduced 
until film boiling started to disappear at the ends. This insured 
that the heat flux was close to minimum. Still pictures of the 
film boiling phenomenon were taken and observations of the temp- 
erature and pressure were also made. The heat flux was again 
increased in steps. At each step a photograph was taken and 
other corresponding observations were made. This procedure 
was continued until the heat flux was such that the wire started 
glowing or the bubble merger mechanism was evident. By 
"bubble merger" we mean that two neighboring bubbles growing 
in the first half of the cycle come into contact with the adjoining 
bubble from the second half of the cycle. Again the heat flux 
was decreased in steps and the process repeated. 
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This procedure was repeated for wires of different sizes 
and various pressures. Each time a new wire was used, the liq- 
uid in the capsule was also replaced. In some cases high speed 
movies were made, to facilitate study of the growth rate of the 
disturbance, with a Hycam motion picture camera at a speed of 
500 frames per second and an exposure time of 1/1250 second. 

3. Data reduction. 

Photographic information was used to make three different 
kinds of observations. We describe them under three separate 
subtitles: 

Wavelength measurements. — Still pictures were printed 
with a 2.54 cm marker as a reference dimension; the wavelengths 
at various heat fluxes were measured. Care was taken to avoid 
situations where bubble merger or deformation of the wave pat- 
tern had taken place. This deformation may be caused by drift 
in the phase angle along the wire or by oscillations generated 
through the movement of the outgoing bubbles. The probable 
error in the measurements of the shortest wavelengths was only 
about ± 4.5 percent. The error was less for larger wavelengths. 

Vapor blanket thickness measurements. — In film boiling, 
a vapor blanket of finite thickness always surrounds the heater 
and there is no liquid contact with the surface of the heater. The 
instability of the upper boundary of this vapor blanket causes a 
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steady release of bubbles. For the application of the theory de- 
veloped in the previous chapter we must obtain the corrected ra- 
dius of the heater wire by adding the minimum blanket thickness 
to R. 

To observe the vapor blanket thickness, representative pic- 
tures of film boiling were enlarged. The minimum diameter of 
the vapor blanket surrounding the wire was measured, using the 
2.54 cm marker again as a reference. This was done to avoid 
any refraction effects that might creep in if the wire diameter 
were chosen as a reference. Furthermore, the wire diameter 
could not be seen clearly through the vapor blanket. The wire 
diameter was subtracted from this measurement to give twice 
the vapor blanket thickness. In other words, we anticipated that 
vapor surrounded the wire equally on the top and the bottom. 
Sometimes a cusp appeared on the blanket in the wake of a de- 
parting bubble, and we have ignored it. 

Measurements of growth rate of disturbance. — Hycam 
movies were viewed on a microfilm viewer. At first, a prelim- 
inary survey was made to sort out those bubbles which grew un- 
disturbed by other adjoining bubbles. Starting with a frame in 
which a bubble had just broken away from the interface, the 
height of the interface was measured from the lowest boundary of 


the vapor blanket. Later, the minimum diameter of the vapor 
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blanket tube surrounding the wire was subtracted from each of 
the above observations to give the amplitude of the wave. 

The minimum height of the interface was used as a refer- 
ence dimension for obtaining dimensionless amplitude. It may 
become more clear if one recalls equation (29) for the amplitude 
of the wave. At a particular value of x, say x=0, one may 
write 

r\ = r\ Q exp(at) (58) 

or 


or 


“T = exp[Q(1.612w t)] 

*L d F 


In (-^) = 0(1.612 ui t) 

TV d cr 


(59) 


(60) 


where is the minimum depth, d^, of vapor blanket and is 

the "most susceptible" frequency for the disturbance in the invis- 

cid fluid in the absence of any geometrical effects of the heater. 

In all cases the dimensionless amplitude, r\/r\ ,was plotted against 

o 

dimensionless time, (1,612 t), on sami logarithmic graph paper. 

F 

The slope of the curve at any instant gave the dimensionless fre- 
quency, 0 . The probable error in the observation of the linear 
growth rate was ± 10 percent. 

Now we would like to compare the experimental observations 


of wavelength and frequency with the theoretical predictions. Be- 
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fore we do this it is necessary to present some sort of correla- 
tion for vapor blanket thickness, because the wire radius has to 
be corrected for it. 


B, Discussion of Results. 

1 . Vapor Blanket Thickness Correlation. 

Baumeister and Hamill [19], while analyzing heat transfer 
from wires in film boiling, developed an expression for vapor 
blanket thickness. Their theoretical model for the film boiling 
configuration was a fairly approximate one; a sequence of spher- 
ical domes connected by annular passages. Neglecting inertia, 
they solved the equations of motion and energy with the assump- 
tion that heat transfer rate is maximum. 

For the heat transfer coefficient for film boiling on wires 
they got: 
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w 



h f9 9(P f~V P 
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(61) 


where a=d^/R. All the physical properties of the vapor in the 

above equation are evaluated at the mean film temperature, i.e. , 

at T=T + AT/2. The latent heat of vaporization has been cor- 
sat 
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re cted to take into account any sensible heat required to heat the 
vapor above its saturation temperature. Assuming a parabolic 
velocity profile and mean film temperature for vapor in the va- 
por blanket, Breen and Westwater [26] gave a corrected expres- 
sion for h : 

fg 




1 + 


0.34 c AT 
v 



] 


Baumeister and Hamill changed the constant in equation (61) 
from 0.23 to 0.35 for small wires (R' ^ 0.05) to take into account 
any heat transfer from the domes. This constant was further 
corrected to 0.485 for larger wires (R' ^ 0.05) to fit the experi- 
mental data. Essentially, the constant in equation (61) was fixed 
from experimental results. This was necessary owing to their 
necessarily simplified model of the film boiling process. 

Their theoretical expression for the vapor blanket thickness 


was 


fg g 
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, 1 / 4 . 


'1 +9R ,2 (1 +A) 2 


] }' 


(62) 


where the term in the second square bracket is close to unity 
for small values of R' and A. 

If heat flux, q, is taken to be an independent variable; the 
temperature difference, AT, between the heater wall and the sur- 


rounding liquid can be written as 
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Using this expression for AT and equation (61) in equation (62) 
gives 


A= exp 
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The group p^/h*^p^ * s a weak function of temperature. We took 
AT=425K, which is typical of the film boiling situation, and 
plotted the present data for cyclohexanol, as well as some old 
data of Sun for acetone [21] against equation (63). The predictions 
were too high and did not correlate the data well for various di- 
mensionless radii. This suggests that we should take another 
look at the various factors influencing vapor blanket thickness. 

The vapor blanket thickness will primarily depend on seven 
independent variables, q, p , h* , g(p -p ), a, p and R. The 

g f g f 3 g 

heat flux q is imagined to be transferred by conduction and used 
fully in the phase transformation. The eight variables, includ- 
ing vapor blanket thickness, d^, can be written in four dimen- 
sions: joule, kg, meter, and sec. Thus, in accordance with 

Buckingham Pi Theorem, we expect the problem to be reducible 
to a relation among four dimensionless groups. For these we 
choose: 


a) The dimensionless vapor blanket thickness, A. 



-41- 


b) The Bond number, Bo, or ratio of buoyant to surface 
tension forces. 

c) A ratio of viscous to surface tension forces. 


B 


qu. 


P K ° 

g fg 


(64) 


d) A ratio of inertial to viscous forces — a kind of Rey- 
nolds number based on the vapor velocity. 


Re = 


qRA 

P h* 

g fg 


(65) 


Baumeister and Hamill's equation (63) involves the first 
three groups but the fourth group is missing because they neg- 
lected inertial terms in their equation of motion. Their assump- 
tion is realistic as long as either the wire radius, or the heat 
flux, or the vapor blanket thickness is small; but the ratio of 
inertial to viscous forces may be of the same order of magnitude 
when the heat flux is high or the vapor blanket thickness is large. 

Radiation losses can also form a significant portion of the 
total heat transfer from a wire in film boiling. We employed 
nichrome wires in all the experiments performed. The total hem- 
ispherical emittance for slightly oxidized nichrome at a tempera- 
ture of 750K has been given by Touloukian [27] as 0.65. Thus 
the radiative heat flux may be written as: 
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4 4 

q = 0.65 o (T - T ) (66) 

r w c 

For AT ~ 425K, which is quite representative for our own vapor 
blanket thickness, equation (66) gives 

5 2 

q^ = 0.132x10 watts/meter (67) 

This heat flux may form as much as 10-20 percent of the total 
heat flux, and is absorbed by the liquid. To know the depth of 
penetration of the radiant heat flux, absorptive properties of li- 
quid cyclohexanol are needed. However, this information is not 
available in the existing literature. Water [28] is a good absorb- 
er of infrared radiation. Values of the absorption or extinction 
coefficient for water as listed in reference [29] indicate that all 
the radiant energy is absorbed in a very thin layer of liquid. 

We presume that this is probably true for cyclohexanol also. The 
radiant energy is absorbed by the liquid adjoining the interface 
and is utilized in generating additional amounts of vapor. 

Depending upon the tocal amount of vapor generated, the va- 
por blanket thickness and vapor velocity in the vapor blanket, 
will adjust themselves in such a way as to optimize the process. 
The Reynolds number, Re, as defined in equation (65) gives a 
relation between vapor velocity and the vapor blanket thickness, 
hence it is an important correlation parameter. 


In Figure 8, we have plotted 21 data points for vapor blank- 
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Correlation for vapor blanket thickness. 
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et thickness on wires in the range 0. 1 1 sr»<,o. 65 during the film 
boiling of acetone and cyclohexanol . While the abscissa is the 
same as Baumeister and Hamill's, the vapor blanket thickness 
has been correlated with the Reynolds number as an additional 
parameter. This correction has been obtained by trial and error 
procedure and as yet, no theoretical reason is offered for taking 
the exponent in the inertia correction term to be 2/3. It may 
also be mentioned that in forming the vapor blanket thickness 
correlation, all of the vapor properties have been evaluated at 
the mean film temperature. 

In all the observations the group, qRis/Ugh*^, varied approx- 
imately from 0.2 to 2.5. Smaller values of this group usually 
occur near the left hand side of the figure. The data scatter 
represents minimum and maximum measurements of vapor blanket 
thickness. The accuracy of these measurements is ±10 percent. 
The data seem to be correlated well by the solid line whose 


governing equation is: 

_1/3 
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( 68 ) 


The constant, 4.60, in equation (68) would correspond to a value 
of 0.32 instead of 0.23 fcr the constant in front of equation (61). 
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2. Comparison of Theoretical Predictions of the 
Wavelength with Experiment. 

Observed wavelengths for cyclohexa.iol boil'ng at a tempera- 
ture of 302. 5K, corresponding to absolute pressure of C.296kPa, 
are displayed in Figure 9. The value of the liquid viscosity 
parameter, M, is close to 5. The theoretical prediction for the 
viscous case is shown in a solid line. The dotted line marks 
the inviscid predictions. The wire radius has been corrected to 
to take into account the thickness of the vapor blanket sur- 
rounding it. 

Although wavelengths were measured for various neat fluxes, 
those displayed in Figure 9 correspond to the lowest heat flux. 
This is done to avoid the longer wavelengths which may be fav- 
ored at higher heat fluxes. We will discuss this point later in 
the text. The data show a wide variability, but the lowest points 
in the range of the data scatter do embrace the theoretical pre- 
dictions which assume that the dominant wavelength is the "most 
susceptible" one. The inviscid predictions would have suggested 
wavelengths about 28 percent too short. 

Figure 10 shows photographs of film boiling for two of the 
data points in Figure 9. Although there is slight phase difference 
across the length of the wire, the Taylor wave is very well de- 
veloped. The bubble release pattern is good. The process is 
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Wavelengths on horizontal cylinders. 
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slower and viscosity seems to have dampened the small interfacial 
disturbances which may be observed in less viscous liquids. 

Wavelength data for M=16 are plotted in Figure 11. Here 
the absolute pressure is 1.06 kPa and the saturation temperature 
is 329K. In this case, observed wavelengths are higher and the 
data again show wide variability. This variability in data could 
be explained from the dispersion relation shown in Figure 5. 

Near the maximum frequency there is a wide region of near- 
neutral stability. Thus, for frequencies slightly tess than the 
maximum, a large range of wavelengths close to X^ is possible. 

3. Wave Growth Rate Analysis. 

Figures 12, 13 and 14 display plots of dimensionless wave 
amplitude versus time on semi-logarithmic coordinates in an in- 
creasing order of R' . In all cases, the liquid viscosity para- 
meter, M, for cyclohexanol is 16. The data contained in these 
figures are for two to three randomly- picked, regularly-growing 
waves during 5 to 10 seconds of motion pictures of film boiling. 

All the figures reveal that bubbles grow linearly during the 
first 12 percent or so of growth. This is the period during which 
we would expect our theoretical predictions of the frequency to be 
valid. Occasionally during the early growth of a bubble, an os- 


cillation of the interface was observed. This was probably caused 
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Fig. 12 Growth of waves on a 0.406mm dia 

wire heater in cyclohexanol. P = l.06kPa, 

q= 0.85 xlO 5 watt/meter 2 , f b = 23 bubble/ 
sec., ai d(r = 45.5 hertz, R' = O.II, Rc=O.I85. 



dimensionless amplitude. 
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Fig. 13. Growth of waves on a 0.51mm dia wire 
heater in cyclohexanol. P= 1.06k Pa, 
q = 0.86 xlO 5 watt/meter 2 , f b =20 
bubble /sec, w dp = 45.5 hertz.R'sO.^.R^O^ 


dimensionless amplitude, t)/tj 
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Fig. 14 Growth of waves on a 0.825 mm dia 

wire heater in cyclohexanol. P = 1 .06 kPa, 
q = I.OIx|0 5 watt /meter 2 , f b = l8 bubble/ 
sec, cu dj _ =45.5 hertz,R'=0.22, R^ = 0.33 
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by the superposition of a disturbance evident in the wake of a 
departing bubble. Such a disturbance is evident in the wake of 
the left-hand bubble in the bottom picture in Figure 10. Data 
corresponding to such disturbances have been identified with dotted 
symbols in the wave growth rate diagrams. One may also note 
that the bubble grows in height to about 75 or 80 percent of the 
wavelength before it leaves the interface. 

Having measured the frequency corresponding to a linear 
growth rate and the dominant wavelength, we next wish to see 
how well the two compare with the predicted dispersion relation. 
In Figures 15, 16 and 17 we trace the dispersion relations for 
the three dimensionless corrected radii and M=16, and we display 
the experimental points on them. The relation between wave- 
length and frequency is borne out quite well in each case. It is 
clear that wavelengths with frequencies slightly less than the 
"most susceptible" frequency can easily occur. In Figure 18, 
observed frequencies are plotted as a function of R ' along with 
the viscous and inviscid predictions. Th^ data fall slightly below 
the viscous predictions, but the inviscid theory would have pre- 
dicted still higher frequencies. 

Wave growth rates for M equal to 52 and 68 are plotted in 
Figures 19 and 20. In both cases, R ' is 0.41. The wavelength 
and frequency are again consistent with the dispersion relation. 
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Fig. 15 Experimental verification of dispersion relation for 

cyclohexanol at 1.06 kPa. q = 0.85xl0 5 watt /meter 2 . 
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Experimental verification of dispersion relation 
for cyclohexanol at 1.06 kPa. 
q = 0.86x10 s watt /meter 2 . 
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Fig. 17 Experimental verification of dispersion relation for 

cyclohexanol at 1.06 k Pa. q = l.0lxl0 5 watt /meter 2 . 
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dimensionless time, (1 612 aj df )t 

F*g. 19 Growth of waves on a 1.03 mm dia 
wire heater in cyclohexanol. P = 4.9 
kPa, q=0.85*l0 5 watt /meter 2 , 
f b = 21 bubble /sec, w dp =46.3 hertz, 
R , = 0.29, Rj, = 0.4 1 . 




dimensionless amplitude, 
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Fig. 20 Growth of waves on a 1.03 mm dia wire 
heater in cyclohexanol. P= 8.05 k pa,q = 

0.98 x I0 5 watt/meter 2 , f b = 22 bubble/sec, 
<^=46.5 hertz, R’ = 0. 295, R^O.41 
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Figures 21 and 22 show bubble growth during the film boil- 
ing of isopropanol on wires of R '=0.172 and 0.535 respectively. 
The liquid viscosity parameter has a value of 264, which means 
that there is very little, if any, effect of viscosity. The data 
presented in these figures have been reduced from high speed 
movies taken by Lienhard and Wong [30]. These movies were 
taken at 2500 ± 40 frames per second. The wave grows linearly 
during the first 15 percent of growth for R^'=0.172 and the first 
35 percent for R c ’=0.535. Figure 22 shows an oscillatory behav- 
ior of the interface during early growth, probably due to superim- 
position of a disturbance caused by the preceding departing bubble. 
We have attempted to ignore the disturbance again and to draw 
the line through the average growth curve. In Figure 23, ob- 
served frequencies for these low viscosity liquids are plotted as 

a function of R ' and compared with theoretical predictions. The 
c 

data tend to be low by only a small percent. 

From the previous observations of growth rates, one may 
also note that the range of non-l’near behavior starts earlier on 
smaller waves than it does on larger ones, and this effect seems 
to be enhanced by viscosity. This may be explained if one looks 
at the drag force experienced by a growing bubble. Keshock and 


Siegel [31] gave an approximate expression for the drag force as: 


cr 

' d 
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dimensionless time, (1-612 

Fig. 21 Growth of waves on a 0.406 mm dia 
wire heater in Iso- propanol. P=98.5 k 
pa, q = 1.7 x 10 5 wait /me ter 2 , f b = 31 

bubble/sec, <u d = 50.0 hertz, R! =0.13, 
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Fig. 22 Growth of waves on a 1.295 mm dia 
wire heater in cyclohexanol. P=98.5 
k Pa , q = 0.5 1 4 x I0 5 watt /meter 2 , 
f b = 2l bubble /sec, w dF = 50.0 hertz, 

R* 5 0.425, Rc=0.54. 
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dominant wavelength on cylinders. 
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where D is an equivalent spherical diameter for the bubble and 
C^is an experimental constant. If in our case, we take n to be 
the equivalent diameter of the bubble, equation (67) can be written 


as: 


F 


d 


dri 

u n — 

f ' dt 


(70) 


Since — is constant during linear growth it may be taken to be 
a function only of the dominant wavelength, i.e. , 


dri 

dt 


uu = f (A.) 


(71) 


For a particular value of r), the drag force increases with liquid 
viscosity and uu; but uu is higher for smaller wavelengths. Thus 
during early growth, small bubbles or bubbles growing in viscous 
liquids may experience more drag. Hence, they approach the 
non-linear growth regime sooner. 


4. Bubble Release Frequency. 

While viewing movies for the bubble growth analysis, we 
also noted the average time after which a bubble was released 
from a node. The bubble release frequency, f^, is found by ">ot- 
ing the number of bubbles released per second. In Figure 24 we 
plot the dimensionless bubble frequency as a function of observed 
wavelength for both the viscous and the inviscid cases. A line 
obeying the conditions: f^O at a. = 0 and df^/dA ->0 at X. = for the 
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flat plate, has been faired through the data. 

The wavelength and liquid viscosity seem to strongly influ- 
ence the ratio of the two frequencies. The "chained" line shows 
the value of f /x^ assumed by Zuber in his prediction of the 
minimum heat flux on flat plates. Zuber obtained an average 
value for this number by averaging the growth rate over ampli- 
tude. The dotted line shows the value of f^/u^p as proposed by 
Berenson. He argued that the growth rate should be averaged 
over the time the bubble takes to grow to its full size and 
breaks off from the interface. Such an average requires full his- 
tory of bubble growth. Berenson did not have this information; 


instead he fixed the value of f^/uu^ em P' r ^ ca ^y fnom the known 


experimental values of minimum heat flux on flat plate heaters. 


Our data suggest an asymptotic value of 0.25 for 


f^/(fl .612) when the liquid viscosity is small and this number 

decreases as liquid viscosity is increased. 


5. Effect of High Volume Fluxes on the 
Dominant Wavelength. 

In the previous sections we saw that frequently wavelengths 
longer than the "most susceptible" wavelength occur in real sys- 
tems. These waves obeyed the dispersion relation well, though 
they did not assume exactly the maximum frequency. As most of 
our experiments were conducted at very low pressures, we pre- 



-67 


sume that the longer, slower- moving waves are preferred when 
large volume fluxes are generated. Large volume fluxes may 
also be caused by high heat f uxes. On the other hand, shorter 
waves may be preferred when small volume fluxes occur. The 
film boiling studies of carbon dioxide near its critical point by 
Abadzic, Grigull and Goldstein [32,33] show the presence of 
nearly stationary waves with wavelengths close to the critical 
value. In this case the volume flux is qjite low. 

In Figure 25 we display dimensionless wavelength as a func- 
tion of volume flux for horizontal cylinders. Wavelength data 
have been obtained both at low pressures and high heat fluxes. 

The coordinates have been reduced to dimensionless inviscid flat 
plate values so that the data obtained at different situations can 
be compared. A definite increase in the dominant wavelength 
seems to take place with volume. Photographs of the film boil- 
ing process at different heat fluxes are shown in Figure 26. 

The longer wavelengths at the larger volume flux in the bottom 
picture can easily be distinguished from the shorter wavelength 
in the top picture. 

6. Some Remarks on the Minimum Heat Flux. 

Various previous authors [18], [20] have pointed out that 


many physical variables, such as material of the wire, the end 



- 68 - 



Fig. 25 Effect of volume flux on the dimensionless dominant 
wavelength. 
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mountings, the material and size of the wireholders, and the 
way the end mountings are attached to the wire, can seriously 
affect the observations of minimum heat flux. Data obtained 
under different conditions have been found to differ with each 
other by as much as 100 to 200 percent. With so many variables 
affecting the true location of the minimum heat flux, an elaborate 
experimental program would be needed to provide reproducible 
and meaningful results. Therefore, we did not make any special 
effort toward measuring the minimum heat flux in the present 
work. 

In Chapter II we made theoretical viscous predictions of the 
minimum heat fluxes on both flat plates and cylindrical heaters, 
and concluded that they should increase with viscosity. As we 
saw in Chapter II, the minimum heat flux depends on the product 
of volume of the bubbles, the number of bubbles per unit heater 
area, and the bubble release frequency. Bubble volume varies 
as dominant wavelength cubed; the area of the heater per bubble 
is proportional to wavelength for cylindrical heaters and to wave- 
length squared for flat plate heaters. Furthermore, the bubble 
release frequency depends mostly on the large period of time 
during which a bubble grows non-linearly . This information 
can only be obtained experimentally. 


Our *v:>eri mental observations show that the bubble release 



frequency decreases with liquid viscosity whereas the dominant 
wavelength increases with it. Hence, at this time it is very dif- 
ficult to say quantitatively what the net effect of viscosity is on 
the minimum heat flux. 

Conclusions. 

1 . An increase of wavelength with the liquid viscosity has 
been measured and found to compare well with the theoretical 
predictions. 

2. The wavelength and frequency measurements are faith- 
ful to the. predicted dispersion relation. 

3. A successful correlation for the vapor blanket thickness 
around a cylindrical heater has been established. 

4. The bubble release frequency seems to depend upon both 
the dominant wavelength and the liquid viscosity. 

5. Generation of large volume fluxes may cause wave- 
lengths that are considerably longer than the "most susceptible" 


wavelength to be favored. 



IV. HYDRODYNAMIC PREDICTION OF THE PEAK HEAT 
FLUX IN INVISCID AND VISCOUS LIQUIDS 
ON LARGE FLAT SURFACES 

Zuber and Tribus explained that the transition from nucleate 
to film boiling occured when vapor jets became Helmholtz unsta- 
ble and collapsed, thereby obstructing the flow of liquid toward 
the heater. In obtaining the expression for the critical velocity 
of these jets, Zuber and Tribus employed the well known result? 
for the interfrcial stability of two parallel-flowing streams of 
inviscid liquids. This instability is caused by pressure fluctua- 
tions and is known as Kelvin-Helmholtz instability. They also 
assumed that the dominant disturbance on these jets was one with 
wavelength equal to critical Rayleigh wavelength. 

In the first part of this chapter, we will review Zuber's 
hydrodynamic model of peak pool boiling heat flux on very large 
flat horizontal surfaces and reconstruct it wherever necessary. 

In the latter pa^t, we discuss the instability of gas jets when 
both gas and liquid are taken to be viscous and use these results 
to obtain expressions for the peak heat flux. 

A. Peak Pcol Boiling Heat Flux When Both the Liquid and Gas 
are Assumed Inviscid . 

f car the peak heat flux the rate of evaporation is very 
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high. The bubble generation frequency is so large that the 
vapor is released in the form of jets. If U is the velocity of 
gas in vapor jets and Aj is the combined area of vapor jets 


supported by an area A of the heater, a simple energy balance 

h 

at the heater surface gives: 


q = 0 h U A /A (72) 

9 fg 9 j h 

The gas velocity at which these jets become Helmholtz unstable 
[34] is given by: 

u = TP- F <73) 

in which X is the dominant wavelength of the disturbance that 
triggers instability of the interface between gas and liquid. 

Zuber also reasoned that on a horizontal surface having no 
geometrical features these jets locate themselves on a two dimen- 
sional square grid spaced on the Taylor unstable wavelength. His 
assumed vapor ►'emoval configuration is shown in Figure 27. 

Zuber could not decide whether the length of the dominant unstable 
Taylor wave should equal the critical wavelength. 


\ c = (? f -P ) 


( 74 ) 


or X , the wavelength which has the fastest growth rate, 
d 
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Fig. 27. Zuber's vapor removal configuration 
on an "infinite" flat plate. 
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X d a 2n vSa/g (p f - P g ^ (75) 

However, later studies of film boiling on cylinders [17, '20] and 

our observations in the preceding chapter have conclusively shown 

that the dominant wavelength is close to X the "most susceptible" 

d 

Taylor wavelength. 

Zuber assumed that the diameter of a vapor jet was half of 
the spacing between the jets. Wong [17] substantiated this assump- 
tion by noting that the bubble diameter was half of the spacing 
between adjacent bubbles in film boiling on cylinders. Thus the 
ratio of area of the vapor jets to the area of heater is a constant: 



The wavelength, X , was assumed by Zuber to be the critical 
Rayleigh unstable wavelength (see ,e.g., [34] page 473 or [35]): 


X = 2rr (X or X )/4 (77) 

H d c 

Substituting equations (73) through (75) in equation (72) we get: 


d 

max 


J/4 


— /p !ag(p -p )]h r 1.196 or 0.S09 ] 

24 9 L fg L J 


(78) 


depending on whether the correct wavelength was X or > 

c d 

Instead, Zuber took the constant in the second square bracket to 
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be unity, which is a good mean value, and obtained equation 

(6) for peak heat flux on a flat plate. 

Subsequent successful hyd rodynam i c predictions of the peak 

heat flux on cylinders [12], spheres [36] and on horizontally and 

vertically oriented ribbons [37] have assumed that the Taylor 

wave fixing the location of jets on the heater is also the "most 

susceptible" disturbance in the jets. The Rayleigh wavelength 

2n X^/4 is longer than the "most susceptible" Taylor wavelength 

X . The photographic observations of references [12] and [36] show 
d 

that the vapor jets on large wires and spheres are too blunt to 

have collapsed by virtue of longer Rayleigh wavelength 2 tt\^/4. 

While the longer wave needs some distance along the jet axis to 

develop, the shorter Taylor wavelength, \ , is already developed 

at the onset and is picked up by the jet. Thus if we propose that 

this wavelength, X , is also the dom lant wavelength in vapor jets 

d 

on a flat plate, Zuber's equation (78) reduces to 

(V) =, - ,4 M'Vf9 V 9 < P '-'Vi ^ 

F inviscid 

= 1.14 q 

max 

Z 

Equation (79) predicts 14% higher peak heat flux on an 
"infinite" flat plate than that predicted by Zuber's equation (6). 
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B. Peak Pool Boiling Heat Flux When Both the Liquid and Gas 
are Assumed Viscous 

In this section we develop an expression for the critical 
velocity of the gas at which the jets become unstable in a vis- 
cous liquid and use this in conjunction with aquation (72) to find 
the peak heat flux. 

The assumed configuration of the gas jet and the liquid 
coU:mn, along with the respective velocity profiles, is shown in 
Figure 28. The inviscid model is also included in the figure 
for the sake of comparison. We will carry out a plane two- 
dimensional analysis, however we will take the effect of transverse 
curvaturo of the vapor jet into account. 

1 . Velocity Profile for the Primary Flow in the Gas Jet and 
Liquid Column. 

Let U denote the velocity of the gas and let U be the 
9 f 

velocity of the lic^jid. At the interface, the gas will tend to drag 
the liquid upward but it has to do so against the force of gravity. 
The gravity force is much larger than the interfacial drag in all 
the situations we are going to consider, hence we assume a zero 
velocity at the interface. The gravity force is balanced by drag 
and the reaction against the liquid at the heater surface. This 
suggests a Poiseuille flow for the gas jet and we accordingly 


take a velocity distribution for the gas as: 
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U = U - U y 2 /R 2 (80) 

g gm gm j 

where. 


R = X /4 (81) 

j d 


The liquid velocity profile has to satisfy the following conditions: 


(i) The velocity at the interface is zero 


U = 0 at y = R 

f j 

(ii) The first derivative of velocity vanishes at the middle 
of the liquid column 

(U.) =0 at y = 2R 

f y j 

(iii) The shear stress at the interface from both the gas 
and the liquid side is the same 


u (U ) = m (U ) at y = R 

f r V q a S 


y g gy j 

(iv) Conservation of mass: The net rate of downward flow 


of liquid is consistant with the amount of vapor generated. 

2rrRj 2tt 2R. 

r r 

d 0 


2rr Rj 2 tt 2Rj 

I I P U ydy d e = [ [ p U ydy 

* * g 9 J J f f 

o o M o R. 


A velocity profile satisfying the above conditions and having 
symmetry about x-axis can be written as: 

Vc 


U = 
f 


u {ri2zri_n ii.fjss. Ha. -s± i.1 \ y */ R 2 

gm '■L 54 ji 81 p J i_ 54 „ 81 p J y j 

f f “ f 


L54 fi f 81 p f J J L 54 /i f 81 p f J j J 


(82) 
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The maximum downward liquid velocity occurs at y = 1.28 Rj and 
may be written as: 


U 


frn 


=-U 

gm 


[19,21 da + 

- 54 4f 


3.14 p g 

— »7 


(83) 


Equation (82) for the liquid profile gives a downward flow of 
liquid near the gas jet and a slight upward flow in the middle. 
This suggests a recirculation of liquid in the liquid column. At 
this point we do not wish to go into details of the liquid motion 
apart from saying that near the peak heat flux the vapor genera- 
tion process is very violent and this type of motion of the liquid 
is possible. 


2. Formulation of the Stability Problem. 

To study the stability of the interface, we superimpose the 

perturbation components of velocity and pressure on the mean 

flow of gas. Let u and v be the perturbation components of 

9 9 

velocity in the x and y directions and p be the perturbation 

9 

pressure. If Pg is the mean pressure in tne gas jet at any 
cross-section, we may write 

u = U + u , v = v and p = P + p 

9 9 9 9 g 

From now on for convenience of typing w*_ do not use the sub- 
script g to denote the physical variables of the gas. The con- 
tinuity equation and equation of motion f or the perturbed flow on the 
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gas side can be written as: 


u + v = 0 

X y 

+ Uu+Uv = - - 


+ V (u 


(84) 

(85) 

x y P 

1 

v + Uv = ~ ~ 

X 

XX 

yy 

p. 

+ >*(v 

+ v ) 

(86) 

t x p 

y 

XX 

yy 


where we have neglected the product of any two small quantities 
and have utilized the fact that the mean flow satisfies the Navier- 
Stokes equations independently. Various conditions these equations 
have to satisfy are: 

(i) The v component of velocity vanishes at the axis of 
symmetry, i.e. at y = 0 

v = 0 at y = 0 (87) 

(ii) The ; irst derivative of the u component of velocity 
vanishes at the axis of symmetry 

u y = 0 at y = 0 (88) 

(iii) The continuity of normal stress at the interface 

2 

-p+ 2uv - ar> /2R. - on = -p at y = R (89) 

y J xx f j 

In the above equation ri (x,t) is the displacement of the inter- 
face from the mean position 

(iv) The continuity of tangential stress at the interface 

u(u + v ) = t at y = R (90) 

y x' f j 

(v) The kinematic condition at the interface can be written 
as 

n = v at y = R (91) 

t 3 

Next, all the variables are made dimensionless _>y choosing 
the radius of the jet, R^, and a characteristic velocity, U Q , as 
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units of length and velocity, respectively. The characteristic 

velocity, U Q , may be taken as the maximum downward velocity 

of the liquid. Keeping the same symbol as for the original 

physical quantity, the dimensionless form of these variables is: 

2 

u = u/U q v = v/U q p = p /pU Q 

x = x/R. y = y/R. t = tU /R 

J J o j 

U = U/U ti = X]/R. T = T/puJ 

o J ° 

Substituting these quantities into equations (84) through (90) we 
get: 


U + V = 

x y 

0 

(84a) 

u + Uu + U V = 

t x y 

-p + R Vu + U ) 
x e xx yy 

-1 

-p + R (v + v ) 
y e v xx yy' 

(85a) 

v + Uv 
t X 

(86a) 


The various boundary conditions become: 


V 

= 0 

at 

y = 

0 

(87a) 

u y 

= 0 

at 

y = 

0 

(88a) 

^e~S*V / ^y +< ^ apU o’V 

II 

-p 

at 

y = 

i 

(89a) 

V v x 

= P. T 
e f 

at 

y = 

i 

(90a) 


where 
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3. Introduction of Perturbation Stream Function. 

We assume that the gas-liquid interface has the following 
two dimensional wave form: 
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i(kx-uit) 


= ri o e 


In the dimensionless form this becomes 


(93) 


„ - _ i a(x-ct) 
T1 > T! o e 


(93a) 


where a is the difnensionless wave number,kRj, and c is the di- 
mensionless wave velocity ou/kU Q . The symbol c denotes a com- 
plex number such that the imaginary part represents the growth 
rate of the disturbance. Now we define a perturbation stream 
function 

iafx - ct) 

(94) 


Y ■ ti ' r (y)e 
o 


such that 


u = Y = fVi 

y 

V =-Y = -iaFn 

X 


(95) 

(96) 


Substituting these definitions of r), u and v from equations (93a), 
(95) and (96) in equations (85a) and (86a) we get: 

-iacF ! n + iaUF 1 ^ - iauVn = -p + r” 1 (-F^o? + F* 1 ^) (85b) 

X © 


-1 2 

-acFri + aFr| = - p^ + (ia Ft) - iFr|) 


(86b) 


Roman superscripts denote differentiations with respect to y. 
Elimination of pressure from equation (85b) and (86b) gives 


F IV -2a 2p - n 


F 11 f a F = iaR e ("(U-c)(F I ^-a 2 F)-U II Fj 


( 97 ) 


Equation (97) is a typical Orr_ Sommerfeld equation and the 
boundary conditions for it are: 
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F (0) = 

0 



(87b) 

F H (0) = 

0 



(88b) 

F 111 -^ 1 = 

icxRe jpf/r]-k7(a 2 - 

-1/2)/(pU 2 R 1- 

U J.j 

-iaR e fcF I (1)-U I F(1)'| 





at y = 1 

(89b) 

F !I + a 2 F = 

t p R /r) 
r e 


at y = 1 

(90b) 


The kinematic condition at the interface can be written as 

c - F (1) (91a) 

4. Solution of Orr-Sommerfeld Equation. 

Our next step is to find a solution of equation (97) with boun- 
dary conditions (87b) through (90b) for values of a close to 1 and 
Re < < 1 . This choice for a is made because the ''most 
susceptible" Taylor wave is also the dominant one on the liquid 
vapor interface. The maximum downward liquid velocity near 
the peak heat flux is very small hence Re < < 1. Thus we write 

a = 1-0 where 8 < 1 (98) 

We also assume a series form for the perturbation function, F 

F = F + 0 F (99) 

o 1 

Substituting equations (98) and (99) for a and F in equation 
(97) and in the boundary conditions (97b) through (90b), and col- 
lecting terms of the order of 8°, etc., we get for terms of 
the order of 9°: 

F IV - 2F 11 + F - 0 (100) 

o o o 
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The boundary conditions for equation (100) are: 

F q (0) = 0 (101) 

F o n (0) = 0 (102) 

F UI (1)-3 F* (1) = i (P+S) (103) 

o u 

Pq^C 1 ) +F 0 < ^ = 1 ^ (104) 

For terms of the order of 8 ^ we get: 

F 1 IV -2F 1 II +F 1 = [(U-c)(F o H -F o ) -u” F q ].iR e /9 (105) 

The boundary conditions on equation (105) are: 

F 1 ( 0 ) = o ( 106 ) 

F ,”(<>) =0 ( i07) 

F^O^SF^O) = -i(P + S)-i £cF^( 1 )-HJ I ( 1 )F q ( 1 ) ] C 108 ) 

F 11 ( 1 )+F 1 (1) =0 ( 109 ) 

In the above equations 

P h p R /n (110) 

f e 

S s a R/(2pU 2 R.) (Ill) 

e o j 

T — t R/t) (H2) 

f e 


The solutions for F and F i can be written as follows: 

o 1 

F o * [o.0434 T - 0.2008 i (P+S) ] (e y -e“ y ) 

+ [0.0732 T + 0.0868 i (P+S) 1 y^-fe^) (113) 
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F = {(P+S)+C [ 0 .5319T-0. 1 478 i(P+S) ]+U ! ( 1 )[o. 3279T-0. 2041 i(P+S)] } 
x{o. 2008 i (ey-e^) - 0.0868i y(e y +e“y)j 
+{2(U-c)[o. 0732T+0. 0868K P+S)]-U I! [o . 0434T-0 . 20081 ( P+S)] ] 
x{y 2 /8 (e y -e~y) } iR e /3 
- { u"[o . 0434T - 0.2008 i (P+S) 1 } 

x {e y (y 3 /24 - y 2 /8) + e*^(y 3 /24 + y 2 /^) } iR e /9 (114) 


5. Kinematic Condition at the Interface. 

The above functions F and F , when evaluated at the inter— 

o 1 

face y = 1 + t] ~ 1, become: 

F o (1) = 0.3279T ~ 0.2041 i (P+S) (115) 

F^l) = 0.2041 i {(P+S)+c [0.5319T-0.1478i(P+S)] 

+ U*( 1) H 0.3279T- 0. 2041 i( P+S)] } 

-c{o.0431T+0.0511i(P+S)} iR^/B 

-U H (1){o.0056T-0.0258i(P+S)} i R e /B (116) 


The wave velocity, c, the perturbation pressure, P, and pertur— 
bation shear, T, are complex quantities. Writing them as the 
sum of real and imaginary parts, we have: 



(119) 
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The kinematic condition (91a) at the interface can be 


written as: 


c = F (1) + S F^l) 


(91b) 


Substituting equations (115) through (119) in (91b) and collecting 
real and imaginary parts we have, 

c = 0.3279 T +0.2041 (1-8) P -9 fo.1085(T.C +T c.) 
r r v i L v i r r i' 

+0 . 0302C.P. -0 . 0302c ^(P +S)+0 . 0669uV 1 )T i -0 . 04 1 6U l ( 1 X P p +S)j 

+R Ic.0056 U n (1) T. - 0.0258 U H (1) (P +5)1 
e(_ t r j 

+Rj0.0431(T c +T c.)— 0.051 1 c.P + 0.0511 c (P + S)] (120) 

e>- i r r r 1 i r r J 


c t - 0.3279 T. -0.2041 (1-0)(P^+S) + 8 fo.0652 (T^c^-T^c.) 

+ 0.0813 (c P.+c P ) -10.0669 (7(1) T +0.0416 l> l P 1 
r i i r r U 

r II — 11—3 

- R e [0.0056 U (1) T r + 0.0258 U (1) P j 

~ R I" 0.0431 (T c -T.c.) -0.0511 (c P.+ c.P ) 1 
eL rrii r i i r J 


( 121 ) 


We are interested in the imaginary part of the wave velocity 
as it represents the rate of growth of the disturbance. For a 
neutral disturbance (c^ = 0),if we put 9 = 0.57(i.e.a= or^= X d ) 
and R « 0 in equation (121) we get: 


0.3279 T. - 0.3204 (P + S) = 0 
1 r 

Equation (122) when written in original variables, becomes 

1.0234 In - !!n. = JL 


( 122 ) 
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6, Normal and Tangential Perturbation Stresses of the Liquid. 

The normal and shear stress perturbations on a rigid 
wavy surface (c = 0) have been given by Brooke Benjamin [38] 
and subsequently used by Craik [39] and Nayfeh and Saric [40] 
for the stability analysis of wind generated waves on thin liquid 
films. The stress perturbations for a boundary layer type of 
external flow with zero pressure gradient are given as: 


fr _ .. 2 

- -a o U 

h f fm 


(124) 


fj 4 00 5/3 2 r jr I nr\ y 1/3 

1 '- 188a <W } < ,25 > 

The liquid velocity profile in our case is more diffuse and 
has a negative pressure gradient in the direction of flow. This 
flow situation could be treated independently but that would require 
making more assumptions and simplification than we would like 
to make. Thus, although the above expressions for perturbation 
stresses do not necessarily represent our flow configuration; 
still it is worthwhile to use the above results and compare the 
final expression for the peak heat flux with experimental data. 
Using equations (124'i and (125) in (123), we have: 

1/3 

[l-i.63^ r /[U^IjRj 2 ]} ]=a/(1R.C f ) (126) 


7. Peak Heat ^lux Prediction. 


Substituting for Rj(= X c j/4) and for U^_ 


from equation (83) in 
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equation (126) and dividing throughout by U 

gc 

tion (73) we obtain 


as given by equa- 


C -0 211 ^[ ,9 - 21 V (54 V + 3- 14 °g /(,8 » f > ] 

gc 

X {,-1.01 t/4 g -(P 9 /» f ) 1/2 • V ( V • M A d ) ] 1 /3 } 

= 2p g /(TT 2 p f A d ) ( 1 2? ) 


On the left hand side of equation (127), the second term in the 
first square bracket is much smaller than the first term at 
pressures much less than the critical pressure, and can be 
neglected. Now for the peak heat flux we write an energy 
balance equation similar to equation (72). 


q = P 
max g 

F 


h (U /2). (tt/16) 
fg P 


(128) 


The use of equation (128) in conjunction with equation (79) in 
equation (127) gives 


< q max p > 

(q max p ) 


0.633 V 


inviscid 


m A 1/2 {l -o. ao7 f _V_ (qma >fr\ n viscid \\ 2 




where M and A are the same as defined in Chapter II and V is 
the product of a gas viscosity parameter similar to M and the 
square root of the ratio of gas and lic^id densities. 
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V h 


p 9 ,/g p f 1/2 p 3/4 

Mg-S'^-Pg) 3/4 


(130) 


Equation (129) is plotted in Figure 29 with M and V as indepen- 
dent variable . For a particular fluid, to every value of M, there 
is a unique value of V. We also trace such a curve for cyclo- 
hexanol in Figure 29. The effect of liquid viscosity is to signi- 
ficantly increase the peak heat flux. 


Conclusions. 

1 . An improved prediction for the peak heat flux for 
nearly inviscid liquids on large horizontal heaters has been 
suggested. 

2. A theoretical analysis for the instability of gas jets in 
viscous liquids has been done and an expression for the peak heat 
flux is obtained as a function of liquid viscosity, gas viscosity 
and density ratio of gas and liquid. 

3. Liquid viscosity seems to strongly increase the peak 


heat flux. 
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V. OBSERVATIONS OF THE PEAK HEAT FLUX 
FOR VISCOUS AND INVISCID LIQUIDS ON 
AN "INFINITE" FLAT PLATE 

A program for the experimental determination of the peak 

heat flux on an infinite flat plate, both at earth normal gravity 

and at elevated gravities, was carried out. By an "infinite" flat 

plate, we mean one with dimensions much larger than the "most 

susceptible" wavelength, X , and one with side walls to avoid any 

d 

inflow of liquid from the sides. The flat plate we used accom- 
modated more than three wavelengths for the alcohols, and about 

two for water, at earth normal gravity; while these numbers in- 
1/2 

creased as g at higher gravities. The observations were made 
for both viscous and nearly-inviscid liquids. Reagent grade liquids 
used were: acetone (CH COCH ), benzene (C H„), isopropanol 

O O O o 

(CH CHOH CH ), methanol (CH„OH), distilled water and cyclohexa- 
3 3 ^ 

nol (CH (CH^ CHOH). 

A. Description of the Experiment. 

1 . Apparatus. 

The flat plate heater and its auxiliary apparatus consisted of 
a heater, a heater support, an emergency shut-off mechanism and 
a reflux condenser. We designed and fabricated each of these 
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elements because commercial models which satisfied our require- 
ments were not commonly available. Hence we will describe 
each of them individually. 

(a) Heater — The flat plate heater -u fulfill the fol- 

lowing requirements. 

The maximum available power without overtaxing the 
existing supply line was 6 kilowatts. This meant that a maximum 
current of 51 amps could be drawn at 110 volts. The heater 
surface area and the heating elements had to be designed for this 
power density. 

The heater was to be capable of operating at either earth 
normal gravity or adaptable to the centrifuge facility for studies 
at elevated gravities. The space available between the centrifuge 
arm and the side wall limited the maximum height to less than 
25 cm. 

To act as an "infinite" flat plate the dimensions of the 
heater surface had to be much larger than the "most susceptible" 
Taylor wavelength. This was necessary to allow enough jets on 
the surface and hence to reduce the end effects. 

The heater design also had to provide one-dimensional heat 
flow to the surface, since two dimensional heat flow would compli 
cate the interpretation of results. 

Apart from being a good conductor of heat, the heater 
material had to be one which would not react with the test liquids 
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With these restrictions in mind the heater was made of pure 
copper with a circular bell shaped configuration. The heater has 
a diameter of 6.45 cm at the boiling surface; 10.16 cm at the 
butt end, and a height of 10.8 cm. A pyrex glass jar allowing 
a liquid head of about 6 cm sits on the heater surface and is 
positioned with the help of two stainless steel flanges. Stainless 
steel has a much lower thermal conductivity than copper and hence 
will have a very small cooling effect on the edges of the boiling 
surface. The total height of the heater with the glass jar is about 
19 cm and it has a final weight of 5 kg. Figure 30 shows a 
sectioned view of the heater. 

The lower portion of the heater contains a central fin and 
two circular rings for transferring heat from the resistance heating 
wire to the main body. Each of the circular rings and central 
fin has square threads machined on its surface to support the 
heating wire. To dissipate 6kw of power, 4 meters of 1.62mm 
dia. nichrome wire was used. These dimensions for the wire 
were chosen so as to minimize the heat input density for the 
maximum available surface area. The wire was covered with 
a fiberglass sleeve to provide electrical insulation from 
the main body of the heater. A thick layer of Sauereisen cement 
was applied on the outside surface to give additional support to 
the wire so that it could not leave the grooves during thermal 
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thermocouples numbered 1-5 



Fig. 30 Sectioned view of the flat plate 
heater. 
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expansion. Being a bad conductor of heat, Sauereisen cement 
also helps in reducing heat loss from the surface. 

To monitor the temperature in various parts of the heater body, 
five 1.2mm diameter holes were drilled. The respective posi- 
tions of the holes and their orientation about the heater axis are 
marked in Figure 31 . Two-hole alumina insulator tubes carrying 
iron constantan thermocouples were pressed into these holes. 

Each of the thermocouple wires measured 0.26mm in diameter 
and had a spark welded junction. The top two thermocouples 
are used to determine the heat flux to the boiling surface. A 
third thermocouple is used for a countercheck on the heat flux 
calculations and a fourth thermocouple to show whether steady 
state conditions are reached. The thermocouple in the heater 
base is used to monitor the temperature for an emergency shut- 
off mechanism . 

(b) Heater Support — A support for the heater is necessary 
for its operation on the ground as well as in the centrifuge. 

Apart from being safe, the heater support had to be light in 
weight and able to function conveniently and efficiently in the 
small space available in the centrifuge. 

Figure 32 shows the heater support. It has two straps with 
a flat base welded inbetween. The material of the support is 
cold-rolled steel. The straps are 6.25cm wide and have a 
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Thermocouples have been designated with 
numbers from I to 5. 

Fig. 31. Orientation of thermocouples about 
heater axis. 
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thickness of 0.95 cm. The strap cross-section was designed to 
accommodate the bearing force encountered at 100g and to incor- 
porate a liberal safety factor. To reduce the weight of the 
support, four large holes 3.18 cm in diameter were bored into 
each of the straps. 

A 1.51cm thick plate with a hole to clear tne outer ring of 
heater wire, and a shoulder to support the heater at the butt end, 
was used as a base of the heater support. An asbestos sleeve 
was placed between the butt end of the heater and shoulder in 
the base plate to reduce conductive heat loss from the heater to 
the support. The lower portion of the heater and the resistance 
wire were finally covered with a shroud made of a thin sheet of 
brass to avoid heat convection to the atmosphere. This heat loss 
could have been significant when the heater is operating in the 
centrifuge . 

(c) Emergency Automatic Shut-Off Mechanism — As a 
la^ge amount of heat is being dissipated over a small area, the 
temperature at the butt end will be rather high. The nichrome 
resistance wire selected for heating purposes melts at about 
1200 K, therefore, we limit our operating temperature at the fin 
end to about 900K. While operating the flat plate heater close to 
the critical heat flux, one may overshoot the actual maximum and 
thus partially insulate the boiling surface. During this period. 
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power is continuously being supplied and it may cause a rapid 
build up of temperature in the body. 

To know the approximate time in which the butt end tempera- 
ture will reach a maximum safe temperature, a one -dimensional 
Schmidt plot was drawn. The Schmidt plot was based on a 20% 
excess power supply. The solution indicates that cnly 20 seconds 
are needed to reach the specified maximum temperature. This 
work was done by Eugene Davis, who was working on the project 
as an engineering aide. Details are given in his laboratory 
report [41]* 

Thus for safe operation of the heater, the power input must 
be terminated quickly. To avoid any delay that may be caused by 
manual observation, it is necessary that some automatic arrange- 
ment should be provided for shutting off the power. 

The output voltage of an iron constanton thermocouple is 
amplified and sent to a Silicon Controlled Rectifier (SCR). When 
the voltage arriving at the SCR reaches a certain preset *'alue, 
a make-and-break switch is activated. It allows the additional 
AC power supply to continue through to a plunger solenoid, 
mechanically linked to the main high amperage circuit breaker. 
When the solenoid is energized it causes the plunger to mc-N/e . 
which in turn pulls open the main power circuit. By turning the 
amplifier, the voltage arriving at the SCR can be made to just 
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reanh the breakover value at the moment the thermocouple 
senses a temperature of 900 K. Figure 33 shows a circuit 
diagram of the shut-off mechanism. A photograph is displayed 
in Figure 34. 

(d) Reflux Condenser — Near the peak heat flux, a tre- 
mendous amount of heat is being dissipated generating large 
amounts of vapor. Simple calculations show that 6 cm of liquid 
head may last only for about two minutes near the critical heat 
flux. Thus we need a very effective refluxing system to capture 
the vapor for reuse. 

For stationary flat plate tests, a glass tube 6.35 cm in dia- 
meter carrying a helical 1.2 cm diameter copper tube was used 
as reflux condenser. The glass tube is 1 meter in height. The 
coil spacing of the copper tube is about 3 cm . The glass tube 
is held with rubber stoppers at the two ends which also support 
the copper tube. Tap water was used as the coolant. 

During tests on the centrifuge, two 15 mm diameter copper 
tubes having fins made of tin sheet were used as a reflux con- 
denser. The fin si 2 e is 28 x 28 x 0.35 mm and the spacing in- 
between fins is 5.5mm. The condenser was designed on the 
basis of forced convection cooling at room temperature. The 
two limbs of the condenser are connected together at one end to 
facilitate connection to a vacuum pump. One may also mention 
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Fig. 33. Circuit diagram of automatic shut-off mechanism. 
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here that the design of the reflux condenser also provided an 
opportunity to make a theoretical analysis of condensation under 
variable gravity. This information has been published in refe- 
rence [42], 

(e) Centrifuge Facility — Elevated gravity tests were per- 
formed on the existing centrifuge facility at the Boiling and Phase 
Change Laboratory. Complete details of the centrifuge design 
are given in reference [213. Figure 35 shows the centrifuge 
designed and built to obtain 100 times earth normal gravity at 
a speed of 360 rpm . When the flat plate heater is put on the 
centrifuge, its boiling surface is 77.7 cms from the center of 
rotation . 

Electrical connections to the centrifuge are made through 
slip rings attached to the drive shaft. These include twelve cii — 
cuits for power input, voltage measurement, and for thermo- 
couples. The drive shaft also carries a tube which is connected 
to the vacuum pump through a stationary vacuum seal. 

A strobe light is placed on a bracket attached to the stand 
supporting the centrifuge. It is triggered by a photo-electric 
pick-off which senses light reflected from a shiny metal piece 
attached to the centrifuge arm . Thus each revolution of the arm 
triggers the light only once. This causes the flat plate heater 
to appear stationary so one can view the boiling phenomenon. 
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Figure 36 shows a schematic diagram of the stationary flat 
plate test and Figure 37 shows a photograph of the flat plate 
heater mounted on centrifuge. 

2. Procedure . 

Before each test, the boiling surface was polished with 

220 grit size emery paper to remove any rust or carbon deposits 

formed during the previous observation. The surface was then 

cleaned with soap and warm water and rinsed with the test liquid. 

This was done to ensure that there was no uncleanliness on the 

surface which might affect the wettability and affect q data. 

max 

The glass jar and reflux condenser were positioned and the 
thermocouple and power connections were made. Before making 
the thermocouple connection to the SCR, a check was made to 
ensure that it functioned properly at the appropriate triggering 
voltage . 

The resistance heating wire was energized and the current 
was gradually increased in steps. At each step the current and 
voltage were noted and thermocouple observations were recorded 
and updated until steady state was reached. It usually took 
about 5 minutes to reach steady state after the current was 
increased to the next higher value. 

Various regimes of nucleate boiling ( e.g . the isolated bubble 
regime and the transition from single bubbles to slugs and 





Fig. 36 
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columns) were easily discernable. The transition from nucleate 
to transitional boiling was identified by noting a slowdown of the 
boiling process and a sudden continuous increase in the ther- 
mocouple reading. The thermocouple reading continued to increase 
for a while, even after the power was removed. In most cases, 
the power was cut off after visual observation of the boiling 
transition, rather than by triggering of the automatic shut-off 
mechanism . 

After removal of power one could also see the boiling pro- 
cess passing over from transitional to nucleate boiling. The 
peak heat flux was computed from the observations of the top 
two thermocouples. Observation of the third thermocouple helped 
in making a countercheck on the heat flux calculations. We 
discuss this and provide an error analysis in appendix A. The 

maximum probable error in q was within 7.7%. 

max 

The peak heat flux observations were made for both viscous 
and inviscid liquids. A few of the observations were repeated to 
check for reproductibility of results. Each time a new observa- 
tion was made the above-mentioned procedure for cleaning the 
boiling surface was used and fresh test liquid was employed. 
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B. Discussion of Results. 

1. Peak Heat Flux on "Infinite" Flat Plates for Nearly Inviscid 
Liquids. 

Peak heat flux data from the present investigations are 
plotted in Figure 38. This shows data only for liquids for which 
liquid viscosity parameter, M, is about 400 or greater, hence 
the viscosity effect is negligibly small. Limited available data 
for clean heater surfaces and reagent grade liquids from other 
sources ( e.g . the ethanol data of Cichelli and Bonilla [8], the 
carbon-tetrachloride and n-pentane data of Berenson [9], and the 
distilled water data of Costello et al. [13]) are also included. 

The peak heat flux has been made dimensionless by dividing 
by Zuber's flat plate value as given by equation (6). The abscissa 
represents the size of the heater in terms of the number of jet 
spacings that can be accommodated on it. This scheme of mark- 
ing the abscissa has beet . chosen to isolate any effect of small- 
er heate^ widths on the peak heat flux. 

From Figure 38 we see that the peak heat flux data are very 
well represented by equation (79) which assumes an "infinite" 
plate model. This mode is apparently valid only for heater 
widths that are greater than about three times the "most suscep- 
tible" wavelength, since the distilled water data corresponding to 
L/X close to 2 lie much below the predicted value. 
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The reason that the water data hang so low could be relat- 
ed to the finite size of the plate. If we look back and go 
through the various assumptions, the assumption most likely to 
be violated is that the ratio of area of the vapor jet to the area 
of the heater is a fixed constant (i «e . rr/16). We further explore 
this matter next. 

2. Peak Heat Flux on Finite Flat Plates for Nearly Inviscid 
Lic^iids. 

In Figure 39, we again plot dimensionless peak heat flux 

against dimensionless heater width for a circular flat plate. 

This time the heat flux is based on the number of jets supported 

by the heater. The choice of a particular integral number of jets 

depends upon the size of the heater and its particular shape. A 

square heater, for example, will support a single jet until its 

width reaches almost twice the dominant wavelength. This 

assumes that there must be a clearance of X^/4 between the jet 

and the side wall to permit a free return of the liquid. The jets 

are also assumed to lie a distance X apart on a sc^jare grid. 

d 

It might be possible for such a square heater to accommodate 

2 jets on a diagonal in the very narrow range of width: 1.707 X^ 

to 2 X . However, this would be an unlikely occurrence. More 
d 

likely, the four-jet grid would establish itself on a slightly 
shortened wavelength when the width gets close to 2 X^, so we 
discount the likelihood of a two-jet configuration. 
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vertical side walls. 
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By the same token, a circular heater will accommodate a 

single jet until the diameter reaches 2 From this point up 

to 2.414 \ a two-jet configuration will be established and at 
d 

2.414 X. , four jets can be accommodated on the plate, and so 
d 

forth. 

As the number of wavelengths that can be accommodated on 
the heater becomes large, it becomes increasingly difficult to 
correctly determine the number of jets supported by the heater 
on account of variations in the dominant wavelength. The 
occurance of a wide spectrum of wavelengths close to the region 
of neutral stability is quite common, as we have observed during 
film boiling on cylinders. 

The effect of these variations should be less and less signi- 
ficant as the number of jets becomes large. Experimental data 
should tend more and more to scatter about the ’’infinite” flat 

plate prediction as we have observed in Figure 38 for L/X 

d 

greater than 3. For smaller plates, the number of jets that can 
be accommodated on the heater can significantly affect the peak 
heat flux. Our distilled water data are very well represented by 
a model which assumes only two standing jets on the heater. In 
Figure 39 curves .m presenting the various number of jets that 
we expect to occur have been terminated at heater sizes where 
the peak heat flux predicted by the preceeding or succeeding 
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curve is the same as for the "infinite” case. 

To further substantiate the above observations and to 
exclude the possibility that the distilled water data are low 
because of the slightly different surface wettability character of 
water as argued by Costello [13], an independent study on square 
flat plates of different sizes and with different liquids was made 
by David Riherd [43], Smooth cold-rolled nichrome ribbons, 
2.29mm thick, and cut to particular sizes, were used as 
heaters. The ribbons were supported f’^om below by ebonite 
sheet and held in brass holders. The brass holders were 
connected to power leads. A rectangular chimney made of tin 
sheet provided side walls above the heaters to avoid any inflow 
of liquid from the sides. A glass window in the chimney was 
used to observe the peak heat flux. 

Riherd's data for various alcohols, and the distilled water 
data point of Costello for a 50.8 x 50.8 mm square plate are re- 
produced in Figure 40. The peak heat flux prediction based on 
the number of jets accommodated on the heater is also included. 
The one-jet model predicts the heat flux well for heater widths 
between one and two Taylor wavelengths. The data start devia- 
ting from the one -jet model as the plate width is decreased below 
one wavelength. This is not surprising. The liquid flow config- 
uration must differ considerably from the assumed potential 
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Fig. 40 Peak heat flux on g square finite flat plate 
with vertical side walls. 
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flow model as the liQjid return route is squeezed into a thin 
annulus near the side wall. 

3. Peak Heat Flux on Flat-Plate and Cylindrical Heaters when 
the Boiled Liquid is Viscous. 

The peak heat flux for cyclohexancl on circular flat plates 
and on cylinders was observed both at earth normal gravity and 
at higher gravities. The data have been nondtmensionalized by 
dividing them by the corresponding inviscid values. 

The comparison of data for cylinders with data for flat 

plates is justifiable if R' is near unity, and we use data in the 

range 0.22 R' < 1.24 here. Sun and Lienhard's [12] inviscid 

hydrodynamic model for the peak heat flux on small cylinders 

assumes the gas jets are separated by X . They took the jet 

d 

radius to be (R+5), where 6 is the thickness of the vapor blan- 
ket surrounding the heater near the peak heat flux. While the 
radius of the jet decreases cs (R + 5), so too does X for small 
cylinders. The effect of a decrease in the liquid return velocity 
due to the smaller jet radius is nearly compensated by this 
shorter Thus this model differs from our assumed flat plate 

model only in that the location of jets is on a line rather than on 
a square grid. 

In Figure 41 are displayed peak heat flux data at earth 
normal gravity. The flat plate data for the smaller values of M 


have been corrected for the effect of finite plate size, since the 
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Fig. 41 Effect of liquid viscosity on peak heat flux for both 
flat plates and cylinders in cyclohexanol at earth 
normal gravity. 
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number of wavelengths accommodated is in the neighborhood of only 
3 after a correction for the viscosity effect is made. 

Now we would like to go back to Chapter IV and take a 
second look at our prediction equation (129). This ec^iation was 
obtained with simplified boundary conditions. In the absence of 
an exact analysis, approximate expressions were used for the 
pressure and shear perturbations on a wavy surface. Thus it 
would be fortunate indeed if equation (129), as it is, represents 
the data well. We find, by comparing the experimental data 
(plotted in Figure 41) with the equation (129), that the predicted 
peak heat flux is higher than the experimental data. However, 
when the constant in the numerator is changed to 0.264, the pre- 
diction fits the data as well as any line could. 

The functional dependence of the peak heat flux on M and V 
is very well brought out by equation (129). From the equations 
of motion, Borishanski [4] recognized that the dimensionless peak 
heat flux of viscous liquids should depend on three other dimen- 
sionless groups. These groups were a liquid viscosity parameter 
N, a ratio of the viscosities of the gas and liquid, and a ratio of 
the densities of gas and liquid. However, from the data for 
slightly viscous liquids his suggested correlation equation (3) was 
only dependent on N. If we return to equation (129) and rewrite 


it in the form 
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(q. 


max. 


(q. 


max 


) 1/2 
F inviscid A d 


/P g (0.793 or empirical constant) 

* ''“f {,- 0 . 94 - 2 - 5 Il^EEeTi 

U d M "f > ... j , j 


Finviscid 
(129) 


where = f(M) and M ~ as described in Chapter II, it 

contains all the groups suggested by Borishanski. Furthermore, 
it is much more general in nature than his equation (3). Equation 
(129) has been derived from purely hydrodynamic considerations. 

Next we test equation (129) with high gravity data for cylin- 
drical heaters and one high gravity data point for a flat plate 
heater. This plot of the peak heat flux as a function of M and V 
is shown in Figure 42. The corresponding predicted values of 
heat flux are connected by lines to the experimental points. The 
data seem to be well represented by equation (129); again using 

0.264 for the constant in the numerator. 


Conclusions. 

1 . The peak heat flux on an "infinite" flat plate for nearly 
inviscid liquids has been found to be 14% higher than predicted by 
Zuber. 

2. The peak heat flux on a finite plate may not be represen- 
ted adequately by the "infinite" flat plate model. The actual 

p 

integral number of. standing jets, rather than one jet per ^ , 



flat plate 
horizontal cylinder 
predicted value 
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Effect of liquid and gas viscosity on the peak heat flux at 
elevated gravities. 
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governs the peak heat flux on smaller plates. 

3. The liquid viscosity tends to increase the peak heat 
flux. The data are correlated by equation (129) with the constant 
in numerator changed to 0.264. 



VI. SUMMARY OF RESULTS 


1 . Liquid viscosity has been shown to increase the "most 
susceptible" Taylor wavelength and decrease the rate of 
growth of the disturbance. The results are not only appli- 
cable to film boiling but to any Taylor unstable processes. 

The effect of liquid viscosity becomes very significant when 
a liquid viscosity parameter, M, is less than 5. 

2. The dominant Taylor wave during the film boiling of cyclo- 
hexanol on cylindrical heaters has been observed. Still 
pictures are made to measure the wavelengths and high speed 
motion pictures to study the wave growth rate. Data for 
wavelength and rate -of -growth are consistant with the theore- 
tical dispersion relation. However, waves slightly longer 
than the "most susceptible" wavelength and having a slightly 
slower growth rate are usually preferred. 

3. A correlation for the thickness of the vapor blanket surroun- 
ding a cylindrical heater has been suggested, based on avail- 
able data for cyclohexanol and acetone. The data are correla- 
ted by equation (68). 

4. Large vapor volume fluxes or high heat fluxes may cause the 
dominant wavelength in p ilm boiling to stretch considerably 
from the "most susceptible" one. 
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5. Zuber's prediction for peak heat flux on "infinite" flat 

plates for nearly inviscid liquids has been slightly modified. 
The suggested expression is 


= 1.14 


^ inviscid 


6. The peak heat flux has been observed on "infinite" flat 
plates at both the earth normal gravity and at highen gra- 
vities and the data have been correlated well by the above 
expression. 

7. The effect of finite size of flat plates on the peak heat flux 
is discussed and some of the related controversial points 
resolved. 


8. The stability of a gas jet in a viscous liquid has been studied. 
A prediction for the peak heat flux is obtained as a function 
of a liquid viscosity parameter, M, and a gas viscosity para- 
meter V. The expression for the peak heat flux is given by 
elation (129). 

9. Equation (129) succeeds in describing the peak heat flux data for 
cyclohexanol on cylindrical and flat plate heaters, at both 
earth normal and elevated gravities when the constant in the 
numerator is changed to 0.264. 

10. The effect of liquid viscosity is to increase the peak heat flux 
greatly when the liquid viscosity parameter.M.is less than 200. 



NOMENCLATURE 



B 


Bo 


c 


c 

v 



h 

rg 

h 


h 

w 

k 

k 

g 

k 

n 

K 

L 


area of the heater 

combined cross-sectional area of vapor jets 
escaping from 

ratio of viscous to surface tension force in vapor 
blanket, q|i g /(p g h f * o) 

the Bond number, ratio of buoyant to surface 
tension forces, r 2 (o f - p ) g/o 

speed of a disturbance in a liquid-vapor interface 

specific heat of vapor 

depth of liquid and vapor phases, respectively 

number of bubbles released per second 

acceleration due to the body force acting in a 
system 

latent heat of vaporization 

latent heat of vaporization corrected for sensible 
heat, h [1 + 0.34 c v A T/h fg ] 2 

convective heat transfer coefficient at the surface 
of a flat plate heater 

heat transfer coefficient for film boiling on wires 
wave number 

thermal conductivity of vapor 
thermal conductivity of heater material 
dimensionless wave number, k - IT)] 

characteristic length of a heater 
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L', R' dimensionless characteristic length of a 

heater, L /g(p_ - p )/a , R' = L' based on L = R 
f 9 

M liquid viscosity parameter, defined in equation (44) 

N Borishanski number, defined in equation (4) 


P 


absolute pressure at the surface of a heater 


P 


P 

9 


P 


f» 


P 

g 


dimensionless perturbation pressure, defined in 
equation (110) 

pressure in a gas jet 

perturbation pressure in liquid and gas, respec- 
tively 


p total perturbation pressure 

o 

q heat flux 


q 

max 


max,. 


peak heat flux for a particular heater geometry 
and for an "infinite" flat plate, respectively 


q 

max^ 


Zuber's prediction as defineo in equation (6) for 
the peak heat flux on an "infinite" flat plate 


q . ,q 

mm min r 


minimum heat flux on a cylindrical heater and on 
an "infinite" flat plate, respectively 


R radius of a cylindrical heater 

R c radius of a cylindrical heater corrected for the 

vapor blanket thickness around it 


R' 


c 


a dimensionless corrected 



R 


e 


Reynold's number as defined in equation (65) and 
equation (92) 


S 

T 

T 


radius of gas jet 

a dirnensionless number, defined in equation (111) 
temperature 

dimensionless perturbation shear stress, defined 
in equation (112) 
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T critical temperature of a fluid 

c 

T saturation temperature of a liquid 

s 

T temperature of heater wall 

w 

f time 


U 


f 


primary velocity of the l»c^iid in a liquid column 



maximum primary downward velocity in the 
liquid 


U 

9 

U 

gro 


U 

f 


» 


u 

9 


V 


v 


primary velocity of the gas in a gas jet 

maximum primary velocity cf the gas 

perturbation velocity of liquid and gas, respecti- 
vely, in the x^direction 

dimensionless gas viscosity parameter, defined in 
equation (130) 

vapor volume flux 


V 


v 

9 


perturbation velocity of liquid and gas, respectively, 
in the y-di recti on 


Greek Letters 


a dimensionless wave number as used in equation 

(93a) 

9 a number less than one 

r dimensionless density, (p - p ) / (p + o ) 

f g f 9 

6 thickness of the vapor blanket on a cylindrical 

heater near the peak heat flux 

A dimensionless vapor blanket thickness in film 

boiling, d /R 
9 

AT temperature difference between a heater wall and 

a liquid 
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h 

A 


A 

d 

A 



s 

i 


Q 

n 

d 

uu 



ordinate of an interface measured from its mean 
position 

dimensionless wavelength, (A/2 -) v / g(5 -p )73:r 

f 9 

= i/k/5 

dimensionless "most susceptible" wavelength 
wavelength 

the "critical" and "most susceptible" wavelengths, 
respectively 

Helmholtz unstable wavelength 

viscosity of liquid and vapor, respectively 

kinematic viscosity of liquid and vapor,respectively 

density of liquid and vapor, respectively 

surface tension between a liquid and its vapor 

perturbation shear stress in liquid and vapor, 
respectively 

perturbation potential function as defined in 
equation (20) 

perturbation function as defined in equation (21) 

perturbation stream function as defined in equa- 
tion (94) 

dimensionless frequency, u u /cr/g® (r> — o ) 

\ f 9 

dimensionless "most susceptible" frequency 
angular frequency 

maximum or "most susceptible" angular frequency 

"most susceptible" angular frequency on a flat 
plate when both the licpjid and gas are considered 

inviscid, /g 3 (p - p )/a/l.612 

V f 9 
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Subscripts 

expt 

refers to experimental data 

F 

"infinite” flat plate 

f 

liquid 

g 

gas 

h 

heater 

i 

imaginary 

j 

jet 

m or max 

maximum 

min 

minimum 

r 

real 
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APPENDIX A 


FLAT PLATE PEAK HEAT FLUX CALCULATIONS 
AND ITS ERROR ANALYSIS 

Heat Flux Calculations. 

The heat flux during any instant of the boiling process was 
evaluated by noting the readings of thermocouples 1 and 2; 
whereas the observation of thermocouple 3 provided a counter- 
check. We again mark the locations of thermocouples 1,2 and 3 
below in Fig. 43 for convenience. All distances are marked in 
millimeters. 



Fie- 43. Location of thermocouples used for heat flux 
calculations. 

For steady state heat conduction in the axial direction, 
Fourier's law can be written as: 

q = -k h H Watts/m 2 (A1) 
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Assuming that the cross-sectional area is constant anc the con- 
ductivity of copper near 100°C is given by [44] ; 


k h (T) = 378.9 - 0.052 (T-100) Watts/mK 


(A2) 


the integration of equation (A1) between 2 and 1 gives: 


9* d = (T -T ) (378.9- .052 - 10 o] \ 

1-2 1-2 2 r l i-2 JJ 


(A3) 


or 


'l-B =( V T P{ 


19890 - 2.73 


( t 1+ t 2 )_ 


L 2 


100 


]} 


(A 3a) 


Similarly, the heat flux between locations (2) and (3) can be 
written as 


f rC T o+ T o) 

q 2 . 3 = (T 3 -T 2 , | 29835 - 4.10 |_ 



(A4) 


For all the observations the contribution of the term in square 
brackets was less than 1.5%. 


Error Analysis. 

The total error in finding the peak heat flux consists of 
three different parts and we discuss each of them under a sepa- 
rate subtitle. 

(i) The Natural Convection Heat Loss. 

This heat loss can affect the observations two ways: 

(a) By giving a radial temperature gradient in the 

heater body: To know this we evaluate the Biot num- 
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ber or the ratio of internal to external resistance. 

For stationary tests, the average natural convective heat 
transfer coefficient from the heater* surface may be taken appro- 

p 

ximately to be 11 watts/(m K). The heat transfer coefficient has 
been calculated by assuming a mean temperature along the heater 
wall and neglecting any curvature effect of the heater. New 


Bi = — ~ 11 x 0 « 0322 ? . = 0.0009 

k h 378.9 

During the centrifuge tests, heat transfer from the outer 
surface of heater takes place due to forced convection. Experi- 
mental data for heat transfer from a cylinder when the flow of 
air is normal to its axis have been obtained by R. Hilpert [45]. 
Using his correlation for heat transfer coefficient as given by 
Kays [46], the Biot number for this situation may be written as: 


Bi = 


25.7 x 0.03225 
378.9 


0.0022 


The Biot numoer in both cases is very small, implyirg thereby 
the outer surface of the heater is nearly insulated and there is 
a negligible heat flow in the radial direction. Hence our assump- 
tion of purely axial heat flow is realistic. 

(b) By convective heat loss from the outer surface of 
the heater between the top thermocouple and the boiling face; 
Energy transferred from the main heater body to the stainless 
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steel flange is eventually lost to the atmosphere through natural 

convection. This heat loss will be considerably less than the 

direct heat loss from the copper surface of the heater. Assuming 

that the heat transfer coefficient is no greater than 25.7 watts/ 
o 

m K and the temperature difference between heater surface and 

5 2 

ambient air is 125K, we obtain 0.014 x 10 watts/m for the 
radial heat loss between the top thermocouple and the boiling 
face. This loss is less than 1% of the lowest heat flux observed, 
(ii) The Error Caused by Measurements: 

The heat flux measurements depend on three 
variables: the thermal conductivity of copper, k^, the difference 

between observations of thermocouples 1 and 2, (T -T^), and 
the distance between the two thermocouples, d^_^. The depen- 
dence of q upon these variables is given by: 


q 


< T 2" T 1> 


1-2 


(A3) 


Th • error in the peak heat flux is the combination of the 
errors in these variables. To analyze the combined effect, we 
differentiate equation (A3) and divide by q. Hence, 


+ A < T 2- T P A < d 1-2> 

CW d 1-2 


Aq _ Ak 

q k 


(A5) 
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The probable error in the measurement of q is: 

ia. 

q y U j K C T 2 - Ti ) ' ^ ' 


(A6) 


The error — can arise from the difference between ac- 
k 

tual conductivity of heater block and the one used in calculations. 
This error should not exceed 1%. 

The thermocouple temperatures were read on a potentio- 
meter calibrated to read the temperature directly. As we are 
interested in the difference of the two thermocouple readings, 
any error in the recorder is cancelled out. The only major 
error is the inaccuracy in reading the scale. The maximum 
value of this error is 6% for the smallest temperature difference 
observed. 

The error in the location of the two thermocouples should 
not be greater than 1%. 

(iii) The Error Due to Occurrence of the Peak Heat Flux 

Before the Steady State Thermocouple Observations are 
Made . 

The main power to the heater was increased in 
steps, waiting each time for steady state to be established. In 
some of the observations, the peak heat flux was observed just 
after the power input to the heater was increased. In such a 
situation, it was difficult to make an accurate reading of the 
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thermocouple while the temperature was still rising. The maxi- 
mum error from this procedure is 4-1/2%. 

The probable error of q can be determined by taking 

max 

the root mean square of all the errors involved. 


* M 


max 


max 


2 2 
(convective error) + (measurement error) + 


(error caused by power input in steps)' 


.2 


r 


C 1 ) 2 + [o) 2 + (6) 2 + (I ) 2 J + (4.5) 2 = 7.7% 



APPENDIX B 


PHYSICAL PROPERTIES OF CYCLOHEXANOL 

To compare the viscous theoretical predictions of the 
Taylor waves with the experimental observations, various phy- 
sical properties of cyclohexanol are needed. The peak heat 
flux calculations also make use of some of these properties. 
Vapor phase properties are needed in the formulation of a corre- 
lation for the vapor blanket thickness. Properties of other test 
liquids, e.g . methanol, acetone, benzene, and iso -propanol, have 
been carefully plotted by Sun [47] . We frequently made use of 
this source in this study. 

The liquid density, p^ t surface tension, a, and vapor pres- 
sure o F cyclohexanol are plotted as a function of temperature in 
Figures 44, 45 and 46. These properties have been obtained 
from references [48] and [49]. The Clausius-Clapeyron equation 
was used for interpolation of vapor pressure data. Viscosity of 
the liquid phase was measured experimentally with a Brookfield 
Syncro - electric Viscometer in the temperature range of 296K to 
335K. These measurements differed a little from the reported 
values. Nevertheless, a mean curve through our experimental 
data and the data reported in [49] has been chosen. Figure 47 
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s a function of temperature. 
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urface tension as a function of temperature 



vapor pressure, k Pa 
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Fig. 46. Vapor pressure of cyclohexanol as a 
function of temperature. 



liquid viscosity /i f x 10 Pa. 
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Fig. 47. Viscosity of cyclohexanol as a 
function of temperature. 




-146- 


shows such a plot for the viscosity of cyclohexarol as a function 
of temperature. The liquid viscosity parameter, M, which is 
nearly equal to the square root of the Borishanski Number, N, 
is plotted as a function of temperature in Figure 48. 

The latent heat of vaporization, h , at various tempera- 

fg 

tures has been calculated by using Watson's relation: 



ty at the normal boiling point is obtained from reference [49]. 

We plot the latent heat as a function of temperature in Figure 

49. The vapor density o at the saturation temperature is cal- 

9 

culated by using the ideal gas law: 

Pg* P/CRT) 

The compressibility factor, 2 = P/p^RT, has been ass ned to 
be unity; which is a very good approximation at low pressures, 
figure 50 shows a plot of the vapor density. 

To obtain the vapor viscosity, use has been made of the 
corresponding state method of Stiel and Thodos. The vapor spe- 
cific heat has been obtained by using Rihani and Doraiswamy's 
method. These methods are treated quite well in reference [50]. 
Figures 51 and 52, respectively, show viscosity and specific heat 



liquid viscosity parameter, 
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Fig. 48. Variation of liquid viscosity parameter, M,for 
cyclohexanol with temperature at earth 
normal gravity. 
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Fig. 49. Latent heat of vaporization of cyclohexanol as a function of 
temperature. 





vapor density, p kg/m 
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Fig. 50. Vapor density of cyclohexanol as a 
function of temperature. 
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Fig. 51 Viscosity of cyclohexanol vapor as a function of 
temperature. 
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as a function of temperature. The _^s viscosity parameter, V, 
for cyclohexanol at ear th normal gravity has been plotted as 
a function of temperature in Figure 53. 

The liquid viscosity parameter, M, and the gas viscosity 
parameter, V, have been evaluated at earth normal gravity. 

These values can be corrected to gravity of interest by dividing 

4/ 

by ./gravity of interest/ earth normal gravity. 



gas viscosity parameter, 
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Fig. 53. Variation of gas viscosity parameter^, 
for cyclohexanol with temperature at 
earth normal gravity. 




APPENDIX C 


TABULATION OF DATA 

The data have been divided into three parts. The first 
part contains wavelength measurements. The second, vapor 
blanket thickness, and the third part, peak heat flux observations. 

In all cases, while evaluating pressure, a correction was 
made for the liquid head acting on the heater surface. This 
correction is quite significant at very low pressures and at high 
gravities. 

For the peak heat flux observations in the centrifuge, the 
frictional head loss in the condenser and other tubing between the 
boiling surface and manometer was also taken into account. 
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TABLE 1 

Wavelength Data of Cyclohexanol Tor Horizontal 
Cylindrical Heaters at Earth Normal Gravity 


R 

(mm) 

R* 



qxlO 5 

fw j 

Observed Wavelength 
(minimum -maximum) 

V ? J 
m 


A 

0.4128 

0.216 

0.296 

5 

.4 

0.S3 

10.41-11.93 

0.50-0.58 



, 

r 





0.96 

11.93-13.97 

0.53-0.67 

0.4128 

0.216 

0.296 

5 

.4 

1.16 

11.18-15.75 

0.54-0.76 

0.5144 

0.268 





0.74 

12.70-14.22 

0.61-0.68 









0.81 

12.70-15.24 

0.61-0.73 









0.82 

14.22-16.76 

0.68-0.80 









0.88 

12.70-13.32 

0.61-0.64 









0.95 

13.46-16.00 

0.65-0.77 









1.00 

14.98-17.53 

0.72-0.S4 

0.6540 

0.343 





0.69 

13.97-17.02 

0.67-0.82 









0.75 

12.70-17.78 

| 0.61-0.86 









! 0.82 

14.48-18.54 

6.70-0.89 



1 






0.84 

13.97-18.54 

0.67-0.89 

0.8000 

0.433 





0.50 

15.49-20.57 

j 0.75-0.99 





1 




0.69 

20.06-21 .08 

: 0.96-1.01 

0.4128 

0.2 

>24 

1.06 

16 

1.11 

11.93-13.72 

| 0.59-0.68 

1 

1 







1.43 

13.72-14.48 

| 0.68-0.72 

0.5144 

0.278 





0.67 

13.46-16.76 

0.67-0.83 

I 

r 







0.62 

13.46-16.76 

0.67-0.83 

0.6540 

0.354 



1 

‘ 

0.72 

j 15.24-18.54 

0.75-0.92 


i 
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TABLE 2 

Vapor Blanket Thickness Data for Horizontal 
Cylindrical Heaters at Earth Normal Gravity 


R 

(mm) 

R' 

qxlO - ^ 

r-x-'s 

m 2 ^ 

QM a 4 

m 

Observed Vapor Blanket 
Thickness 

(minimum - maximum) 

. A 1 V 

o h *o 
9 fg 

h* 

9 fg 1 

i 

4 i 

| dg(mm) 

A= dg/R 

Cyclohexanol 

0.2667 

0.144 

0.86 

10 

1.25 

! 0.12-0.15 

0.45-C.55 

0.4128 

0.216 

1.16 

64 

3.42] 0.12-0.16 

0.28-0.38 

1 

| 


1 

1.70 

104 

5.55 j 0.18-0.21 

0.44-0.51 



f 

2.12 

130 

6.86! 0.19-0.23 

0.47-0.56 

1 


0.227 | 

1.01 

11.6 

2.30 

0.17-0.24 

0.40-0.58 

0.5144 

0.268 ; 

0.96 1 

53 

3.60 

0.13-0.17 

0.26-0.34 



\ 

1.68 

92 

6.00 

, 0.21-0.26 

0.41-0.51 



0.278 

0.67 

7.3 

1.90 

0.15-0.19 

0.29-0.37 



0.287 

0.79 

3.4 

2.64 

0.17-0.22 

0.33-0.43 



0.294 

0.79 

1.8 

2.54] 

0.11-0.13 

0.22-0.26 



1 

1.02 

2.4 

3.30] 

0.10-0.13 

0.20-0.26 

0.6540 

0.343 

0.84 

! 46 

3.86 

0. 13—0. 17 

0.20-0.26 


1 

1 

0.92 ] 

! 50 

i 

| 4.25 

i 

0.21-0.26 

0.32*0.40 

Acetone + 

0.2667 

0.161 

1.04 

0.36 | 

1.44] 

0.04-0.05 I 

0.15-0.19 

0.3200 

0.203 

1.30 

0.45 ! 

2.27 

0.04-0.06 

0.13-0.19 

0.4128 

0.257 

1.19 

0.42 

2.63 

0.07-0.08 

0.17-G.18 

0.6440 

0.409 

1.01 

0.35 

3.55 

0.05-0.06 

] 0.08-0.1 

0.6540 

0.458 

0.88 

0.31 

3.10 

0.07-0. 10 

0.11-0.16 

0.8250 

0.514 

1.17 

0.41 

5.16 

0.09-0.11 ; 

0.11-0.14 

1.0280 

0.650 

1.11 

0.39 

6.20, 
1 

0.06-0.09 

i 

0.06-0.09 

1 


+ 


Data reduced from observations of Sun as reported by Lienhard 
and Carter [21] . 
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TABLE 3B 

Peak Heat Flux Data of Cyclohexanol 
for Horizontal Cylindrical Heaters 


R 

(mm) 

| R ’ 

g/g 

e 

P 

(kPa) 

M 

V 

W' 10 '* 

( 5 ) 

Smax 

(<1 

m ox . . . 

inviscid 

0.4128 

0.216 

1 

0.30 

5.4 

125 

3.60 

8.00 

0.5144 

0.268 

1 


0.30 

5.4 

125 

2.72 

j 6.43 

0.6540 

0.343 



! 0.30 

5.4 

125 

2.74 

6.88 

0.8001 

0.419 



0.30 

5.4 

125 

3.31 

8.60 

0.5144 

0.284 



2.68 

29 

312 

4.60 

4.15 


1 j 

0.289 



5.00 

44 

366 

4.89 

3.43 


f 

0.294 



7.72 

66 

| 425 

4.79 

2.88 

0.6E 

>40 

0.354 j 


1 

1.06 

i 16 

233 

3.18 

4.50 



0.362 



2,41 

j 27 

300 

3.59 

3.58 



0.368 



5.47 

47 

375 

4.00 

2.90 



0.373 



7.17 

61 

420 

4.26 

2.83 

0.4128 

0.222 



0.94 

14 

215 

3.25 

4.34 



0.228 



2.76 

30 

317 

3.37 

2.87 



0.233 



7.86 

67 

429 

3.69 

2.01 


i 

0.241 



13.31 

104 

500 

4.29 

1.92 



0.669 

8. 

5 

3.54 

20 

196 

5.27 

3.01 



0.997 

18. 

.3 

5.34 j 

23 

181 

6.24 

2.66 



1.194 

25. 

7 

6.89 

26 

182 

8.67 

2.96 



0.680 

8. 

,3 

8.00 

37 

247 

5.48 

2.23 



1.020 

18. 

,3 

10.41 

37 

218 

7.16 

2.24 



1.220 

25. 

7 

11.99 

36 

209 

8.10 

2.29 



0.672 

8. 

,3 

5.45 

28 

224 

5.71 

2.76 

i 


0.998 

17. 

,8 

8.00 

30 

205 

7.00 

2.50 



1.207 

25. 

,7 

9.44 

32 

193 

7.95 

2.44 



0.687 

8. 

,3 

10.40 

44 

262 

5.80 

2.13 



1.011 

17. 

,8 

11.58 

39 

221 

6.69 

2.08 



1.220 

25. 

.7 

12.82 

39 

209 

8.10 

2.19 



0.702 

8. 

3 

16.52 

63 

295 

5.48 

1.65 



1.032 

17. 

8 

17.99 

55 

251 

7.44 

1.90 

! 


1.240 

25. 

.7 

19.23 

53 

239 

8.10 

1.90 









APPENDIX D 


ON THE USE OF SI UNITS IN THIS STUDY 

The international scientific and technical community is grad- 
ually accepting one common system of measures known as the 
International System of Units (SI) as a substitute to the various 
customary units. Recognizing this trend, we have presented all 
the results in this work in SI units. For the convenience of those 
who may wish to use this information in other units, a table of 
conversion factors from SI units to Metric or English is given 
below. 


Quantity 

SI Units 

Equivalent 



Metric Units 

English Units 

density 

kg/m 3 

, 3 

1 kg/m 

0. 06241 bm/cuft 

surface tension 

Newton/m 

1 kg/s 2 

0.0685 Ibf/ft 

viscosity 

Pa* s 

(N:s/m 2 ) 

I 

f 1 kg/s-m l 

1000.0 centipoises 

1 

! 413.22 Ibm/ft hr 

| 

pressure 

i 

k Pa j 

101.97 kg/m 2 
7.5006 torrs 

I 0.1451 psi 
| 0.01 bars 

energy flux 

W/m 2 

2.39x10 kg-cal/m s 

| 0.3171 Btu/ft 2 hr 

: 

specific heat 

kJAg *K 

i 

0.2389 kg-cal/kg°c 

j 0.2389 Btu/lbm^ 

latent heat 

j 

kj/kg J 

0.2389 kg-cal/kg 

l 0.4303 Btu/lbm 
1 
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